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ABSTRACT 

Simplified  finite  element  modeling  of  a  stiffened  cylinder  subjected  to 
underwater  explosion  was  investigated.  The  use  of  smearing  the  stiffeners  into  the 
base  structure  as  well  as  beam  modeling  using  SOR  (Surface  of  Revolution)  beam 
elements  were  used  in  the  model  simplification  process.  The  dynamic  response  and 
overall  global  deformation  were  then  compared  between  the  fully  discretized  stiffened 
cylinder  model  and  the  simplified  models.  The  study  first  examined  the 
effectiveness  of  smearing  stiffeners  onto  a  flat  plate.  The  smearing  of  stiffeners  into 
a  cylindrical  shell  orthotropically  was  then  examined.  Next,  beam  modeling  of  both 
unstiffened  and  stiffened  cylinders  was  investigated.  Finally  an  integrated  beam/shell 
model  of  a  stiffened  cylinder  was  developed.  These  models  were  subjected  to  the 
same  underwater  explosive  loading  for  numerical  study.  The  analysis  showed  that 
when  comparing  the  dynamic  responses  caused  by  underwater  explosions  between 
the  discrete  model,  the  beam  model,  and  the  beam/shell  model  of  a  stiffened  cylinder, 
globally  similar  results  could  be  produced. 
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I.  INTRODUCTION 

Detailed  finite  element  modeling  of  a  submerged  structure  including  internal 
structures  subjected  to  underwater  explosion  is  costly.  The  cost  can  be  equated  in  both 
computational  time  and  more  importantly  manpower.  As  a  result,  a  simple,  but  still 
accurate  model  of  a  structure  was  investigated.  The  objective  of  the  study  was  to  develop 
an  efficient  means  of  representing  the  stiffened  cylinder  while  still  providing  accurate 
dynamic  response  predictions  to  underwater  explosion.  Past  work  at  the  Naval  Postgraduate 
School  focused  on  comparison  of  numerical  modeling  to  the  physical  testing.  Cylindrical 
models  have  been  previously  subjected  to  various  types  of  underwater  explosions  and  their 
results  have  been  closely  approximated  by  the  computer  code.  Fox,  Kwon,  and  Shin 
[Ref.  1]  and  Chisum  [Ref.  2]  have  all  demonstrated  that  the  finite  element  and  boundary 
element  programs  closely  approximate  experimental  results.  Their  research  work  clearly 
indicates  that  limited  numerical  modeling  can  be  conducted  without  the  need  to  physically 
build  the  actual  structure  for  experimental  testing. 

The  studies  [Ref.  1]  did  not  explore  modeling  cylindrical  shells  with  internal 
structures.  This  type  of  modeling  involves  adding  internal  equipment  and  structures  to  the 
outer  shell,  which  dramatically  increases  the  degrees  of  freedom  of  the  problem.  This 
directly  correlates  to  a  long  computational  time  for  solutions.  This  study  will  examine 
two  means  by  which  a  simplified  model  and  discrete  model  can  be  numerically  tested  to 
produce  similar  dynamic  responses  due  to  an  underwater  explosion.  The  first  technique 
models  the  stiffened  cylinder  with  SOR(Surface  of  Revolution)  beam  elements  and  the 
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second  technique  involves  smearing  the  stiffeners  onto  the  cylindrical  shell  structure. 

Using  a  material  smearing  technique  is  prevalent  in  composites  where  smearing 
material  properties  reduces  the  complexity  of  the  structure  when  subjected  to  various  forms 
of  loading.  Raftenberg  [Ref.  3]  demonstrated  that  smearing  in  structures  required  that  the 
structure  be  entirely  axiymetric  or  have  regions  which  varied  circumferentially  about  its 
axis  of  symmetry  with  rapid  periodicity.  The  smearing  technique  was  used  for  static  loading 
only.  Pitarresi,  Caletka,  Caldwell ,  and  Smith  used  a  similar  finite  element  material-smearing 
technique  to  dynamically  test  a  printed  circuit  board  with  favorable  results  [Ref.  4].  Both 
studies  dramatically  reduced  both  complexity  and  computational  time  with  their 
simplification  methods. 

The  finite  element  and  boundary  element  programs  utilized  for  providing  the 
numerical  analysis  for  this  study  are  called  USA/DYNA3D.  VEC/DYNA3D  is  a  finite 
element  analysis  code  developed  by  Lawrence  Livermore  National  Laboratory  [Ref.  5].  The 
Underwater  Shock  Analysis,  USA  [Ref.  6  ]  is  a  boundary  element  code  that  is  based  on  the 
Doubly  Asymptotic  Approximation  (DAA)  which  was  developed  by  Geers  [Ref.  7].  The 
linkage  between  the  two  codes  was  developed  in  1991  and  has  provided  acceptable  results 
of  dynamic  responses  of  cylindrical  shells. 

One  of  the  applications  of  the  USA  code  developed  by  DeRuntz  [Ref.  6]  allows  the 
beam  to  be  modeled  with  SOR  (Surface  of  Revolution)  elements.  The  SOR  element 
approximates  a  beam  element  with  a  surface  of  revolution  element.  The  SOR  element 
contains  a  number  of  sub-elements  around  the  circumference  of  each  surface  of  revolution 
element.     The  SOR  branch  interpolates   sub-elements  using  a  circumferential  harmonic 


function  which  ultimately  reduces  the  number  of  surface  of  revolution  degrees  of  freedom. 
The  original  discrete  model  of  the  cylinder  analyzed  in  this  study  was  provided  by 
Protection  &  Weapons  Effects  Department,  Code  67  Carderock  Division  of  the  Naval 
Surface  Warfare  Center  and  is  provided  in  Figure  1.  This  cylinder  is  called  the  ONR 
cylinder  and  is  a  test  vehicle  which  is  currently  being  planned  for  a  series  of  underwater 
explosions.  This  study  will  first  examine  the  effects  of  smearing  stiffeners  onto  a  flat  plate. 
Next,  the  discrete  model  will  be  modeled  without  internal  stiffeners  with  Belytschko-Lin- 
Tsay  4  node  shell  elements  due  to  its  computational  efficiency  in  DYNA3D  program.  The 
unstiffened  model  will  be  modeled  as  a  beam  with  SOR  elements  and  an  integrated 
beam(SOR)/shell  model  will  be  developed.  All  three  models  will  be  subjected  to  simulated 
underwater  explosions  and  their  responses  will  be  analyzed  and  discussed. 

The  discrete  model  with  internal  stiffeners  will  then  be  modeled.  Based  on  the  study 
of  smearing  stiffeners  into  a  flat  plate,  the  ONR  stiffeners  are  smeared  into  the  shell  of  the 
cylinder  orthotropically.  Beam  model  with  SOR  elements  will  be  developed  as  well  as 
integrated  beam/shell  model.  The  internal  stiffeners  will  be  incorporated  into  beam  model 
by  adjusting  section  properties  of  beam.  The  discrete  and  two  simplified  models  will  be 
numerically  analyzed  by  applying  the  explosive  charge  first  at  the  symmetric  point  of  all 
three  models.  Detailed  analysis  will  be  provided  at  beam/shell  interface  as  well  as 
endplates  and  deep  frame  stiffeners  to  compare  discrete  model  to  simplified  models.  The 
three  models  will  then  be  exposed  to  an  off-center  charge  and  analysis  will  be  performed  as 
previously  noted. 


Figure  1.  Naval  Warfare  Surface  Center's  ONR  Cylinder 


4 


n.  UNDERWATER  EXPLOSION  CHARACTERISTICS 

The  underwater  explosion  testing  that  is  simulated  by  USA  computer  code  provides 
the  user  many  options.  For  this  study,  however,  the  same  options  and  explosion 
characteristics  were  used  so  that  a  simlilar  comparison  between  the  discrete  or  full  model 
and  simplified  models  could  be  made.  All  models  developed  in  this  study  were  subjected 
to  an  underwater  shock  wave  modeled  as  a  spherical  wave  with  a  peak  amplitude  of  2585 
psi  and  a  decay  constant  of  0.3493  msec.  This  approximation  to  the  actual  shock  wave  uses 
an  exponential  pressure  relationship  and  is  good  for  pressure  greater  than  about  1/3  of  the 
peak  value  [Ref.  8].  This  peak  pressure  pulse  is  generated  using  60  lbs  charge  of  HBX-1 
explosive  at  a  10  foot  depth  and  a  standoff  distance  of  27.2  feet.  The  attack  geometry  is 
provided  in  Figure  2.  The  pressure  profile  is  provided  in  Figure  3  for  a  reference. 

The  10  foot  underwater  depth  was  selected  in  order  to  fall  within  the  gas  bubble 
radius  of  15.8  feet.  This  allows  the  gas  bubble  to  vent  to  the  surface  and  not  migrate  to  the 
surface  of  the  cylindrical  surface.  If  the  bubble  were  not  to  vent  to  the  surface  then  a 
significant  amount  of  explosive  energy  would  be  subjected  to  the  structure(approximately 
47%)  [Ref.  8].  This  was  done  in  order  to  eliminate  the  effects  of  bubble  pulse  on  the 
simulations  perfomed  in  this  study.     Free  surface  effects  were  also  ignored  in  this  study. 


HBX-1  charge 


Figure  2.  Attack  Geometry 


60#  HBX-1  UNDEX  Pressure  Profile 


3000 


0.5  1 

time  in  (msec) 


Figure  3.  Pressure  Profile  for  HBX-1  60  Pound  Charge 


m.  FLAT  PLATE  SMEARING 

The  first  area  of  model  simplification  that  was  studied  involved  examining  the  effects 
of  smearing  stiffeners  onto  a  flat  plate.  The  aim  of  the  study  was  to  find  an  optimal  way  to 
smear  the  stiffeners  without  affecting  the  global  dynamic  response  of  the  flat  plate.  From 
this  study,  the  objective  was  to  use  these  results  for  finding  an  efficient  means  of  smearing 
the  stiffeners  onto  the  cylindrical  shell. 

The  plates  used  in  this  study  were  either  clamped  or  simply  supported  at  the  edges 
and  were  subjected  to  a  unit  step  pressure  load  of  lOOOpsi.  All  plates  examined  were  88" 
x  88"  with  a  stiffener  spacing  of  4  inches.  The  stiffeners  are  1.25"  thick  and  0.25"  wide  and 
were  modeled  in  both  directions.  A  comparison  between  the  discrete  and  smeared  models 
was  accomplished  using  VEC/DYNA3D  finite  element  programs.  The  discrete  model  had 
the  stiffeners  fully  modeled  to  the  base  plate  whereas  the  smeared  model  had  the  stiffeners 
smeared  onto  the  base  plate.  The  smeared  value  was  simply  an  adjustment  in  the  material 
property  of  the  base  material.  The  adjustment  is  mad  within  the  material  definition  portion 
of  INGRID.  The  smeared  value  was  obtained  using  a  well  known  formula  derived  from 
Timoshenko  and  Woinowsky-Kreiger  [Ref.  9]: 


12  EI   (1-v2) 

E'  -  E  +  ^ 

ah3 


where 

E'=  fictitious  smeared  elastic  modulus 
E  =  elastic  modulus  of  plate  and  stiffener 
Is  =  moment  of  inertia  of  the  stiffener 
a  =  spacing  between  two  stiffeners 
h  =  thickness  of  the  plate 
v  =  Poison's  ratio 

The  density  of  the  base  material  was  also  adjusted  to  take  into  account  the  mass  of  the 

stiffeners. 

The  discrete  model  of  a  0. 1875"  thick  plate  with  clamped  edges  was  generated  and 
subjected  to  the  unit  step  pressure  loading.  Figure  4  is  the  dynamic  response  of  the  discrete 
model  at  0.5  msec.  The  transverse  deformation  plot  is  magnified  by  a  scale  factor  of  10  for 
clearer  visualization  of  the  deformation.  Next  the  smeared  value  was  computed  for  a  0. 1875 
"  thick  plate  with  the  given  parameters  for  the  stiffeners.  The  smeared  model  was  subjected 
to  the  same  pressure  loading  and  Figure  5  shows  the  deformation  of  the  smeared  model.  The 
fictitious  smearing  value  was  approximately  85  times  larger  than  the  original  material's 
elastic  modulus  ( EV  85*E). 

This  modeling  shows  the  case  where  the  stiffeners  dominated  the  overall  response 
of  the  original  structure.  Notice  that  when  comparing  the  two  results,  it  was  observed  that 
the  discrete  model  had  a  much  larger  deformation  gradient  near  edges  of  the  plate.  The  area 
near  the  edges  of  the  discrete  model  governs  the  global  deformation  of  the  structure  and  for 
discussion  purposes  will  be  called  the  critically  deformed  area.  The  glaring  difference  in 
global  deformation  between  the  discrete  and  smeared  models  can  be  accounted  for  by  the  low 
stiffness  of  the  discrete  model  base  plate.  In  other  words,  the  original  discrete  model  has 
1/85  of  the  globally  smeared  model's  stiffness  in  the  base  plate.  As  a  result,  stiffeners  in  the 
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Figure  4.  Global  Deformation  Plot  of  Discrete  Model  of  0.1875"  Plate  With  Clamped 
Edges,  Magnified  by  a  factor  of  10,  taken  at  0.5msec 


Figure  5.  Globally  Smeared  Model  of  0.1875"  Plate  With  Clamped  Edges 
Magnified  by  a  factor  of  10,  taken  at  0.5msec 


mid-section  of  the  discrete  model's  base  plate  are  less  constrained  from  the  side  edge 
constraints.  This  produces  a  very  large  gradient  near  the  edges  of  the  base  plate  and  a 
minimum  gradient  at  the  mid-section  of  the  plate.  The  contrapositive  of  this  analysis  is  that 
the  smeared  plate  is  uniformly  constrained  in  both  directions  therefore  the  globally  smeared 
plate  has  a  uniform  deformation  gradient  over  the  entire  plate.  Although  in  Figures  4,5  the 
deformation  is  shown  at  0.5  msec,  the  comparison  is  very  similar  even  at  3.0  msec. 

The  most  effective  way  of  simplifying  the  discrete  model  and  obtaining  the  proper 
deformation  gradient  at  the  critically  deformed  area  is  to  smear  locally.  That  is,  leave  the 
discrete  stiffeners  modeled  around  the  edges  of  the  plate  in  the  critically  deformed  area, 
however  locally  smear  the  remaining  stiffeners  into  the  base  plate  as  in  the  globally  smeared 
model  using  the  same  smearing  value  as  in  the  globally  smeared  model.  Figure  6  shows  the 
dynamic  response  of  the  locally  smeared  versus  the  discrete  model.  The  local  smearing  did 
not  affect  the  critically  deformed  area.  As  a  result  the  dynamic  response  of  the  locally 
smeared  plate  was  in  excellent  agreement  with  the  discrete  model. 

A  second  study  of  the  same  stiffened  plate  was  undertaken  using  simply  supported 
boundary  conditions  at  the  plate  edges.  The  dynamic  responses  produced  similar  results  to 
what  was  observed  for  the  clamped  plate  modeling.  The  simply  supported  0.1875"  thick 
plate  models  (discrete  and  globally  smeared)  are  provided  in  Figure  7.  The  results  were 
very  similar  out  to  3.0  msec  so  only  0.5  msec  is  shown. 

Next,  a  study  was  conducted  to  analyze  the  case  where  the  ratio  of  fictitious  elastic 
modulus  to  original  elastic  modulus  is  lowered.  By  increasing  the  plate  thickness  to  1",  the 
ratio  was  reduced  to  2  (E^  2E).  This  time,  the  results     were  expected  to  show  excellent 
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6a)  Locally  Smeared  Model 


6b)  Discrete  Model 

Figure  6a,b.    Comparison  of  Locally  Smeared  model  of  0.1875"  Plate  versus  Discrete 
Model  of  0.1875"  Plate.  Both  With  Clamped  Edges  and  with 
magnification  factor  of  10,  taken  at  0.5  msec 
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(a)  Discrete  Model  with  Simply  Supported  Edges 


(b)  Globally  Smeared  Model 

Figure  7.  0.1875"  Plates  With  Simply  Supported  Edges:  (a)  Discrete  Model  (b) 

Globally  Smeared  Model.  Both  with  10X  magnification  ,  at  0.5  msec 
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agreement  between  the  discrete  and  globally  smeared  models  due  to  the  increased  stiffness 
in  the  base  plate.  Again  the  discrete  model  was  generated  and  subjected  to  the  same  unit 
step  pressure  loading.  Figure  8  shows  the  result  of  the  1"  discrete  plate  model  versus  the 
globally  smeared  model  with  both  models  having  clamped  edges.  Figure  9  compares  both 
discrete  and  globally  smeared  1"  plates  with  simply  supported  edges.  As  predicted  the 
responses  were  very  close  to  each  other  for  both  sets  of  boundary  conditions. 

Based  on  the  study  conducted  on  smearing  flat  plates,  the  smearing  technique 
directly  supported  the  well  known  and  established  formula  derived  in  Ref.  9.  The  present 
study  provided  an  additional  criteria  for  determining  the  proper  approach  to  smearing 
stiffeners  onto  a  base  structure  based  on  the  original  plate's  thickness  and  elastic  modulus. 
If  the  ratio  between  fictitious  elastic  modulus  and  original  elastic  modulus  is  too  high 
erroneous  solutions  will  be  obtained  if  global  smearing  is  undertaken.  If  the  ratio  is  small, 
then  global  smearing  is  an  effective  means  for  simplifying  the  structure. 
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(a)  1"  Discrete  Model  With  Clamped  Edges 


(b)  Globally  Smeared  1"  Plate  Model 

Figure  8.  1"  Plates  With  Clamped  Edges:  (a)  Discrete  Model  (b)  Globally  Smeared 
Model.   Both  with  10X  magnification  factor  and  taken  at  0.5msec. 
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(a)  Discrete  Model  With  Simply  Supported  Edges 


(b)  Globally  Smeared  Model 

Figure  9.    1"  Plates  With  simply  Supported  Edges:  (a)  Discrete  Model  (b)Globally 
Smeared  Model.  Both  at  10X  magnification  factor  and  at  time  0.5  msec 
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IV.  MODELING  OF  UNSTIFFENED  ONR  CYLINDER 

When  first  considering  the  task  of  simplifying  a  stiffened  cylinder,  a  detailed  study 
was  undertaken  to  investigate  simplification  techniques  in  general.  Two  techniques  were 
investigated  for  modeling  internal  structures  such  as  stiffeners  into  the  external  shell.  The 
technique  previously  investigated  studied  smearing  the  stiffeners  onto  the  base  structure. 
This  knowledge  will  be  used  when  investigating  orthotropically  smearing  the  stiffeners  into 
the  cylindrical  shell  later  in  this  report.  The  second  simplification  technique  involves  using 
SOR  elements  to  effectively  model  the  cylindrical  shell  and  internal  structures.  This  section 
of  the  work  details  how  to  model  an  unstiffened  cylinder  with  SOR  elements.  The  SOR 
beam  model  will  be  developed  and  compared  to  the  unstiffened  shell.  Finally,  an  integrated 
SOR  beam/shell  model  of  the  unstiffened  cylinder  will  be  generated  and  all  three  models  will 
be  subjected  to  the  same  explosive  loading. 
A.         DISCRETE  UNSTIFFENED  ONR  CYLINDER  MODEL 

The  three-dimensional  model  is  shown  in  Figure  10.  Its  dimensions  conform  to 
those  depicted  on  Figure  1  except  for  no  internal  stiffeners.  This  particular  model  uses  HY- 
100  steel  for  the  shell  and  mild  steel  for  the  endplates.  Both  materials  for  modeling 
purposes  were  considered  an  isotropic  elastic  material.  Material  properties  used  for  this 
model  are  provided  in  Table  1 . 
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Mild  Steel 

HY-100 

Density 

(lbm/ft3) 

490.0 

490.0 

Poisson's  Ratio 

0.3 

0.3 

Yield  Strength 

(ksi) 

32.0 

108.0 

Young's  Modulus 

(psi) 

2.9x1 07 

2.9xl07 

Tangential  Hardening 

Modulus 

(psi) 

5.1xl04 

5.02xl04 

Table  1.  Material  Properties  for  ONR  Cylinder 


x 

A 


y   <• 


Figure  10.  ONR  Cylinder  Modeled  Without  Internal  Stiffeners 
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Although  the  model  shown  is  a  full  model,  a  half  model  of  the  cylinder  was  used  in 
the  numerical  analysis  to  take  advantage  of  symmetry  condition.  The  yz-plane  was  the 
symmetric  plane  from  Figure  10.  Even  with  a  half  cylinder  model,  there  are  20  elements 
around  the  half  circumference  and  45  elements  in  the  axial  direction  with  an  element  length 
of  5.625  inches.  The  element  length  represents  the  spacing  between  the  shallow  frame 
stiffeners. 

The  discrete  model  was  subjected  to  an  underwater  shock  with  USA/DYNA3D 
programs  with  shot  geometry  shown  on  Figure  2.  The  major  response  of  the  cylinder  to  the 
shock  is  in  the  yz-plane.  Since  there  was  a  boundary  condition  imposed  on  yz-plane,  there 
is  no  global  motion  in  the  x-axis.  As  a  result,  only  global  deformation  responses  in  the  z- 
direction  are  compared  in  Figure  1 1  at  discrete  states  of  time.  The  plots  are  magnified  by 
a  factor  of  50  for  improved  visual  ease.  At  early  time,  the  deformation  appears  symmetric 
except  for  slight  pinching  near  endplates.  The  endplates  behave  similarly  to  the  clamped 
edges  on  the  flat  plates  in  that  a  large  deformation  gradient  forms  around  each  endplate. 

Another  way  to  analyze  the  dynamic  response  is  to  plot  the  displacement,  velocity 
and  acceleration  at  desired  points  on  the  cylinder.  A  comparison  of  the  axial  displacement, 
velocity  and  acceleration  will  indicate  the  cylinder's  accordion  motion.  Plots  in  the  z- 
direction  indicate  the  whipping  motion  of  the  cylinder.  These  time  history  plots  will  be 
presented  later  in  this  chapter  to  compare  the  discrete  model's  response  to  the  simplified 
model's  dynamic  response. 
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11a)  State  time  of  0.5  msec 


lib)  State  time  of  1.0  msec 


lie)  State  time  of  1.5  msec 


Figure  lla,b,c.  Deformation  Results  of  Unstiffened  ONR  Cylinder  Model 
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lid)  State  time  of  2.0  msec 


lie)  State  time  of  2.5  msec 


llf)  State  time  of  3.0  msec 


Figure  lld,e,f.    Deformation  Results  For  Unstiffened  ONR  Cylinder  Model  (cont.) 
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B.         SOR  BEAM  MODEL 

The  SOR  beam  element  formulation  is  based  on  approximating  the  shell's  curvature 
with  surfaces  of  revolution.  This  type  of  approximation  has  extensive  applications  for 
modeling  containers,  tanks,  domes  and  has  been  done  for  a  long  period  of  time  [Ref.  9]. 
A  small  displacement  of  a  point  on  the  surface  of  revolution  can  be  resolved  into  two 
components.  One  of  the  components  is  in  the  direction  tangent  to  the  curvature  of  the 
surface  about  its  SOR  axis.  The  other  component  is  in  the  direction  of  the  normal  to  the 
beam  surface.  The  displacement  about  the  surfaces  of  revolution  are  estimated  by  applying 
a  sine  and  cosine  harmonic  function  for  approximating  the  response  on  the  surfaces  of 
revolution.  The  harmonic  function  is  then  incorporated  into  the  double  asymptotic 
approximation  procedure  for  estimating  the  fluid  structure  interaction  in  an  infinite  medium. 
As  previously  noted,  DeRuntz  [Ref.  6]  developed  an  application  to  USA  codes  that  assigns 
a  harmonic  function  to  each  of  the  SOR  elements  or  branches.  The  harmonic  function  will 
describe  the  SOR  fluid  degrees  of  freedom  and  more  closely  predict  the  responses  to  the 
underwater  explosive  loading.  During  this  study,  the  beam  was  modeled  with  separate  sine 
and  cosine  functions,  however  it  was  proven  that  a  more  accurate  numerical  result  was 
obtained  when  both  sine  and  cosine  functions  were  used  for  modeling  the  beam  with  SOR 
elements.  The  net  result  was  that  the  unstiffened  and  stiffened  cylinder  could  be  effectively 
modeled  as  a  beam  using  the  SOR  element  technique. 

The  SOR  beam  model  was  modeled  with  the  same  number  of  finite  elements  (45)  that 
was  used  to  model  the  discrete  cylinder  along  its  axis.  The  model  is  shown  in  Figure  12  and 
was  subjected  to  explosive  loading.     The  beam  had  the  same  material  properties  as  the 
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cylinder's  shell  and  its  cross-sectional  area  and  moments  of  inertia  were  computed  using  the 
cylinder's  dimensions.    That  is,  the  cross-section  area  of  the  SOR  beam  is 

A  -  ^di  -  4) 

4 
where  d0  and  dt  are  the  outside  and  inside  diameters,  respectively.   Similarly  the  second 
moment  of  inertia  and  polar  moment  of  inertia  of  the  cross-section  are 


64 


and 


j  ■  £tf  -  *t) 
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Figure  12.  SOR  Beam  Model  With  45  Elements 
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The  global  deformation  of  the  beam  model  was  plotted  for  state  times  0.5  msec  to  3.0 
msec  and  is  shown  in  Figure  13.  The  plots  are  magnified  by  a  factor  of  50  for  enhanced 
visual  recognition.  There  is  an  appearance  of  some  crumpling  or  wrinkling  at  a  few  peaks 
with  an  increase  in  time.  The  deformation  gradient  at  the  ends  are  not  as  distinct  as  with  the 
cylinder,  however  the  overall  shape  and  curvature  is  in  excellent  agreement  with  the  discrete 
model.  The  displacement,  velocity,  and  acceleration  time  histories  were  also  generated  and 
will  be  discussed  later  in  the  chapter. 


24 


13a)  State  time  of  0.5  msec 


-i   i   ' 


13b)  State  time  of  1.0  msec 


1  i  i  '    ~ 


13c)  State  time  of  1.5  msec 


13d)  State  time  of  2.0  msec 


3     I  ^SE 


13e)  State  time  of  2.5  msec 


1        I        '        '        I        I        '        ■        '        ■        '        I        I        I  I        l        r 


13f)  State  time  of  3.0  msec 


Figure  13a,b,c,d,e,f.     Deformation  Plots  for  SOR  Beam  Model  with 

50X  magnification  factor 
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C.        INTEGRATED  SOR  BEAM/SHELL  MODEL 

A  third  study  was  undertaken  to  simplify  the  original  unstiffened  cylinder  by 
integrating  SOR  beam  elements  with  shell  elements.  This  type  of  model  would  not  only 
reduce  the  size  of  the  discrete  model  but  would  now  be  sufficiently  small  enough  to  allow 
for  additional  internal  structures  to  be  added  to  the  shell  section  of  the  beam/shell  model. 
The  same  harmonic  functions  for  the  SOR  elements  were  used  and  the  shell  elements  were 
generated  in  the  same  manner  as  in  the  discrete  model.  The  success  of  this  model  solely 
depended  on  developing  the  proper  interface  between  SOR  beam  and  shell  elements  that 
would  produce  accurate  dynamic  response  to  the  explosions.  The  first  consideration  for  the 
interface  between  both  types  of  elements  was  using  internal  beams  to  connect  SOR  elements 
to  shell  elements.  This  consideration  proved  to  be  totally  ineffective  because  there  was 
entirely  too  much  local  deformation  at  the  beam/shell  interface. 

The  next  type  of  interface  explored  was  the  use  of  a  dry  endplate  with  thickness  equal 
to  the  thickness  of  the  shell.  More  favorable  dynamic  responses  were  obtained,  however 
there  were  sizeable  differences  in  the  axial  motions  of  the  model.  This  type  of  model  is 
represented  in  Figure  14.  The  model  is  shown  as  a  full  model,  even  though  a  half  model  was 
used  to  reduce  computational  time  and  take  advantage  of  symmetry  conditions. 

1.         Variation  of  Material  Properties  for  Beam/Shell  Interface 

A  sensitivity  study  was  undertaken  by  changing  the  material  property  of  the  dry 
endplate  being  used  as  the  interface.  Both  the  thickness  of  the  plate  and  Young's  modulus 
were  adjusted  in  order  to  observe  changes  in  the  deformation,  both  axially  and  in  the 
transverse  or  z-direction.  The  transverse  or  whipping  motion  of  the  model  was  not  affected 
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Figure  14.        Integrated  SOR  Beam/Shell  Model  with  Dry  Endplate  Used  as 
Interface 


by  large  changes  in  the  material  properties  for  the  endplate.  The  axially  direction,  however 
was  affected  by  changes  in  either  endplate  thickness  or  Young's  modulus. 

It  was  noted  that  the  dynamic  responses  were  more  sensitive  to  changes  in  endplate 
thickness.  Only  when  Young's  modulus  was  adjusted  in  a  large  amount,  there  was  a 
noticeable  change  in  the  dynamic  response  of  the  entire  model.    After  performing  numerous 
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numerical  studies  using  the  dry  endplate,  it  was  determined  that  the  dry  endplate  produced 
favorable  results  for  the  modeling  the  integrated  SOR  beam/shell  model.  Furthermore,  the 
endplate's  material  properties  can  be  approximated  by  the  following  formula: 


3*  P*  a4        x 

<,  o 


4*  E*  t3 


where: 


P=  max  incident  pressure(psi) 

a=  radius  of  plate 

E=  Young's  modulus 

t=  thickness  of  plate 

6  =  tolerance  of  axial  displacement  for  the  endplate 

The  formula  was  simplified  from  the  exact  solution  for  the  deflection  of  a  circular  plate 

loaded  at  the  center  [Ref.  9].  Results  from  application  of  the  formula  have  produced  an 

integrated  model  that  has  very  similar  whipping  motion  responses  and  slight  differences  in 

axial  direction  responses.  The  reasons  for  slight  differences  in  axial  direction  responses  can 

be  attributed  to  approximating  the  shell  with  SOR  elements.    Although  SOR  elements 

produce  good  whipping  motions,  the  axial  motions  are  slightly  off  because  the  shell  elements 

not  only  bend  axially  but  also  bend  in  the  radial  direction  causing  differences  in  the  overall 

result.  The  SOR  beam  model  did  not  adjust  locally  for  the  mass  of  the  endplates  which  also 

causes  differences  in  the  axial  direction.  Locally  smearing  the  properties  of  the  endplates 

will  be  examined  in  the  discussion  on  stiffened  beam/shell  models. 

Figure  15  shows  a  comparison  of   the  SOR/shell  model  with  a  3/16"  thick  dry 

endplate.  Figure  15a  shows  the  model's  global  deformation  at  0.5  msec  after  the  shock  has 
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hit  the  structure  and  before  the  interface  endplate's  material  properties  were  adjusted. 
Figure  15b  shows  the  global  deformation  after  the  endplate's  material  property  was  adjusted. 
Notice  that  the  unadjusted  model  has  a  noticeable  amount  of  axial  deformation  at  the 
endplate  as  compared  to  the  adjusted  model  which  shows  very  little  axial  deformation  at 
endplate.  The  tolerance  for  the  adjusted  model  was  set  at  .01  inch.  Using  a  small  tolerance 
ensures  that  there  will  be  very  little  deflection  of  the  dry  endplate  which  improves  the  overall 
global  deformation  of  the  model.  Also,  the  small  tolerance  ensures  a  smooth  transition  of 
deformation  at  the  interface.  There  is  some  pinching  of  the  shell  caused  by  the  rigidity  of  the 
dry  endplate,  however  it  will  be  seen  later  in  this  chapter  that  the  overall  global  deformation 
and  dynamic  responses  at  the  beam/shell  interface  are  in  good  agreement  with  the  discrete 
model.  This  type  of  model  also  provides  for  a  smooth  transition  of  the  dynamic  responses 
from  beam  to  shell  at  the  interface  region. 
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15a)  3/16"  Thick  Endplate-  No  Adjustment  to  Material  Properties 
50  X  Magnification  Factor  ,  state  time  of  0.5  msec 


15b)  3/16"  Thick  Endplate,  Material  Property  Adjusted, 
50  X  Magnification  Factor,  state  time  0.5  msec 


Figure  15a,b.  Global  Deformation  Plots  for  SOR  Beam/Shell  Model: 
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The  beam/shell  model  with  adjusted  material  properties  provides  very  similar 
deformation  results  when  compared  to  the  discrete  model,  however  there  is  one  disadvantage 
to  this  type  of  modeling.  Adjusting  the  elastic  modulus  causes  a  longer  computation  time 
due  to  the  increased  stiffness  in  each  element.  The  time  step  size  used  in  the  explicit  finite 
element  computer  program  is  computed  by  the  following  formula: 


At 


SL 


\ 


E 
P 


where: 


S=  safety  factor  (use  number  less  than  1) 
L=length  of  smallest  element  in  any  direction 
E=  Young's  modulus 
p=  material  density 

By  selecting  a  larger  Young's  modulus,  the  overall  time  step  is  decreased  which  increased 

the  overall  computational  time. 

After  formulating  the  proper  way  to  integrate  SOR  beam/shell  elements,  a  study  was 

done  to  find  a  dry  endplate  that  would  produce  accurate  results  but  reduce  the  computational 

time.    A  thicker  endplate  was  first  examined  since  it  would  help  reduce  the  ratio  of  Young's 

modulus  to  material  density.      The  Young's  modulus  for  the  1/2"  thick  endplate  was 

calculated  using  the  endplate  formula.    A  1/2"  thick  dry  endplate  was  selected  for  the 

integrated  SOR/shell  model  since  this  produced  an  adjusted  Young'  modulus  of  about  100 

times  less  than  the  3/16"  thick  endplate's  Young's  modulus.    Both  models  were  subjected  to 

the  same  explosive  testing  with  the  charge  location   at  the  symmetric  point  on  the  shell. 
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Figures  16-18  compare  the  axial  displacements,  velocities  and  accelerations  between  the 
3/16"  and  1/2"  thick  endplate  models  out  to  3.0  msec.  Since  only  elastic  response  was 
examined,  the  responses  out  to  3.0  msec  produces  the  most  effective  time  to  compare  the 
results.  There  is  exceptional  agreement  with  all  three  plots,  however  there  is  a  slight 
difference  in  magnitude  and  frequency  which  can  be  accounted  for  by  different  masses  of  the 
endplates.  The  1/2"  thick  plate  produced  slightly  larger  displacements  in  the  axially 
direction.  Figures  19-21  compares  the  z-direction  or  whipping  motion  responses  of  the  two 
models  at  the  shock  point  on  both  models.  There  is  very  tight  agreement  between  both 
models  for  the  displacement,  velocity,  and  acceleration.  This  clearly  suggests  that  the 
estimation  of  the  material  property  of  the  endplates  is  a  very  effective  tool  for  integrating  the 
SOR  Beam/Shell  model.  Even  though  a  thicker  endplate  was  used  the  overall  time  step 
was  (9.027xl0'8)  which  is  30  times  less  than  the  original  time  step  of  (2.5X10"6).  Further 
optimization  of  this  problem  would  produce  a  dry  endplate  of  about  1.75"  that  would  have 
a  time  step  of  about  80%  of  the  original  time  step.  Because  DYNA3D  code  allows  a  material 
choice  option  for  modeling  a  material  as  a  rigid  body,  the  time  step  problem  can  be 
eliminated.    This  material  option  will  now  be  discussed. 
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Axial-Displacement  For  Two  Beam/Shell  Adjusted  Models 
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Figure  16.  Axial  Displacement  of  Two  Integrated  SOR  Beam/Shell  Models 
Both  With  Material  Property  Adjustment  to  Endplates 


33 


Axial-Velocity  For  Two  Beam/Shell  Adjusted  Models 
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Figure  17.  Axial  Velocity  Response  for  Two  Integrated  SOR  Beam/Shell  Models 
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Axial-Acceleration  For  Two  Beam/Shell  Adjusted  Models 
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Figure  18.  Axial  Acceleration  Response  for  Two  Integrated  SOR  Beam/Shell  Models 
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Z-Displacement  For  Two  Beam/Shell  Adjusted  Models 
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Figure  19.  Whipping  Displacement  Response  for  Two  Integrated  Models 
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Z- Velocity  For  Two  Beam  Shell  Adjusted  Models 
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Figure  20.  Whipping  Velocity  Response  for  Two  Integrated  Models 
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Figure  21.  Whipping  Acceleration  Response  for  Two  Integrated  Models 
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The  rigid  body  material  option  on  DYNA3D  is  built  into  the  computer  code  and 
only  requires  three  material  property  inputs  (E,  p,  v  ).  The  thickness  of  the  endplate  was 
first  selected  at  3/16"  and  then  varied  to  see  how  well  this  material  option  would  predict  the 
dynamic  responses.  It  is  recommended  that  the  thickness  of  the  endplate  be  very  close  to 
the  thickness  of  the  cylindrical  shell  or  too  much  local  deformation  will  occur  at  the  interface 
area  due  to  the  increased  mass  distribution  of  the  thicker  endplate.  The  structure  created 
with  this  material  option  will  then  behave  as  a  rigid  body  dynamically  when  subjected  to  any 
type  of  loading.  The  SOR  beam/shell  model  was  redesigned  using  this  material  option  and 
tested.  Figure  22  shows  the  global  deformation  plots  for  the  model  using  the  rigid  body 
material  option.  As  with  the  previous  modeling  of  the  beam/shell  interface,  the  rigid  body 
model  produces  pinching  at  the  endplate  interface.  Figure  23  shows  a  comparison  plot 
between  the  adjusted  3/16"  thick  endplate  model  and  the  beam/shell  model  with  the  rigid 
body  endplate  at  0.5  msec.  Both  models  were  in  excellent  agreement  with  each  other.  The 
big  advantage  of  using  this  material  option  is  that  the  original  values  of  Young's  modulus  and 
density  can  be  used  which  in  return  does  not  affect  the  original  time  step.  This  significantly 
reduces  the  overall  computational  time  of  the  model  and  makes  it  a  viable  choice  for 
simplifying  the  cylindrical  shell. 
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Figure  22a)  State  time  of  0.5  msec 


Figure  22b)  State  time  of  1.0  msec 


Figure  22c)  State  time  of  1.5  msec 

Figure  22a,b,c.  Deformation  Results  for  Integrated  SOR  Beam/Shell  Model  Using 
Rigid  Body  Material  Option  with  DYNA3D 


38 


Figure  22d)  State  time  of  2.0  msec 


Figure  22e)  State  time  of  2.5  msec 


Figure  22f)  State  time  of  3.0  msec 
Figure  22d,e,f.  Deformation  Results  for  Integrated  SOR  Beam/Shell  Model  (cont.) 
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Figure  23a)    Deformation  Plot  of  3/16"  Thick  Endplate  for  SOR  Beam/shell  Model 


Figure  23b)  Deformation  Results  of  Rigid  Body  Endplate  for  SOR  Beam/shell  Model 


Figure  23  a,b.  Deformation  Results  for  3/16"  Thick  Endplate  SOR  Beam/shell  Model 
Versus  Rigid  Body  Endplate  SOR  Beam/shell  Model  of  Unstiffened 
ONR  Cylinder 
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D.         COMPARISON  OF  DISCRETE  TO  SIMPLIFIED  MODELS 

This  section  will  compare  the  results  of  the  explosive  loading  on  the  discrete 
unstiffened  model  to  those  of  the  simplified  SOR  beam  and  integrated  SOR  beam/shell 
models.  Comparisons  will  be  done  by  plotting  the  displacement,  velocity  and  acceleration 
in  both  the  axial  and  trans verse(z)  directions  for  all  three  models  at  selected  points.  The  z- 
direction  or  whipping  motion  response  will  be  examined  at  the  center  point  nearest  to  the 
charge.  The  axial  or  accordion  motion  responses  will  be  monitored  at  the  left  endplate  for 
all  three  models. 

For  the  cylinder,  the  center  point  nearest  to  the  charge  (so-called  side  shell  node)  and 
at  the  top  of  the  cylinder  at  the  center  (called  top  shell  node)  and  center  point  away  from 
charge  (called  back  side  node)  could  be  selected  for  comparison.  Because  the  cylinder  has 
both  global  and  local  deformations  at  the  three  nodal  points,  the  locations  are  expected  to 
have  different  responses.  Depending  on  which  motion  is  dominant,  the  difference  may  be 
large  or  small.  If  the  global  response  is  dominant  over  the  local  response,  the  three  points 
will  have  a  similar  response.  Otherwise,  the  three  points  will  have  quite  different  responses. 
In  the  z-direction,  the  global  response  dominanted  over  the  local  response.  This  can  be 
clearly  seen  in  Figure  24  which  compares  the  three  shell  nodal  points  to  the  same  beam  nodal 
point  on  the  simplified  SOR  beam  model.  Notice  how  the  responses  are  very  similar  and  that 
the  top  and  back  side  shell  nodes  have  a  slight  delay  before  the  deformation  occurs.  This  is 
due  to  the  shock  wave  progressing  through  the  cylinder  where  it  hits  the  front  side  shell  node 
first.  Except  for  the  delay,  the  responses  are  quite  similar  in  shape  and  magnitude.  This 
study  will  use  the  front  side  shell  node  only  for  comparison  purposes. 
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Figure  24.  Whipping  Displacement  Plot  Showing  Comparison  of  Shell  Node  Response 
for  Unstiffened  Shell  versus  SOR  Beam  Model 
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The  axial  displacement,  velocity,  and  acceleration  time  history  plots  for  the  discrete 
and  two  simplified  models  are  provided  in  Figures  25-27.  The  beam/shell  model  was 
modeled  with  the  rigid  body  material  option  using  a  3/16"  thick  endplate  for  more  accurate 
results  and  a  quicker  computational  time.  There  is  however  noticeable  differences  between 
all  three  models  for  all  three  responses.  An  explanation  for  the  differences  is  that  the  SOR 
beam  model  will  not  have  the  same  local  nodal  point  response  as  the  same  point  on  the 
cylinder.  The  local  response  in  the  axial  direction  is  more  dominant  than  the  global  response. 
The  beam/shell  model  and  the  beam  model  did  not  account  locally  for  the  mass  of  the 
endplate  which  causes  some  error  in  the  axial  direction.  If  the  magnitude  of  the  local 
responses  are  examined,  there  is  about  a  45%  difference  between  shell  and  beam  model 
displacements.  Figure  28  shows  the  global  deformation  plot  at  0.5  msec  for  all  three  models. 
The  overall  global  deformation  is  in  good  agreement  except  for  local  deformation  at  the 
beam/shell  interface. 

Figures  29-3 1  provide  the  whipping  responses  for  all  three  models  at  the  location  of 
the  closest  point  to  the  charge.  These  plots  show  excellent  agreement  between  all  three 
models  with  the  error  less  than  10%.  The  SOR  approximation  predicts  responses  better  in 
the  transverse  direction  than  in  the  axial  direction. 
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Figure  25.  Axial  Displacement  for  All  Three  Models 
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Axial-Velocity  for  All  Three  ONR  Models 
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Figure  26.  Axial  Velocity  for  All  Three  Models 


Axial-Acceleration  for  All  Three  ONR  Models 
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Figure  27.  Axial  Acceleration  Plot  for  All  Three  Models 
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28a)  Global  Deformation  for  Unstiffened  Cylinder 


2Sb)  Global  Deformation  for  SOR  Beam  Model 


28c)  Global  Deformation  Result  For  SOR  Beam/Shell  Model 


Figure  28a,b,c.  Global  Deformation  Plots  for  All  Three  Models  at  0.5msec 
with  50X  Magnification  Factor 
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Figure  29.  Whipping  Displacement  for  All  Three  Models 
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Figure  30.  Whipping  Velocity  Plot  for  All  Three  Models 
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Figure  31.  Whipping  Accelearation  Plot  for  All  Three  Models 
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And  finally,  Table  2  is  provided  to  compare  the  size  of  the  model  with  the  overall 
computational  time  for  each  model.  The  SOR  beam/shell  model  data  is  for  the  rigid  body 
material  option  for  modeling  the  beam/shell  interface.  Clearly,  one  must  make  a  decision 
as  to  whether  to  use  SOR  beam  elements  or  original  shell  elements.  Clearly  the  accuracy  lies 
in  using  shell  elements,  however  if  a  section  of  the  cylinder  model  requires  isolation  for 
detailed  internal  modeling,  then  the  beam/shell  model  can  be  used  to  simplify  the  design  and 
to  reduce  computational  time.  The  next  part  of  the  study  will  examine  modeling  the  ONR 
stiffnened  cylinder. 


Unstiffened  ONR 
Cylinder 

SOR  Beam  Model 

Integrated 
Beam/shell  Model 

Number  of 
Nodes 

1140 

47 

534 

Number  of  Wet 
Elements 

1080 

45 

360 

Computational 
Time  (sec) 

1968 

6.44 

742 

Table  2.  Comparison  Data  For  all  Three  Models  of  Unstiffened  ONR  Cylinder 
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V.  MODELING  OF  STIFFENED  ONR  CYLINDER 

In  this  section  the  simplification  of  the  ONR  cylinder  is  investigated.  The  original 
ONR  cylinder  with  stiffeners  will  be  modeled  first.  Next,  the  stiffener  smearing  technique 
will  be  used  to  model  the  cylinder  with  orthotropic  shell  elements.  The  SOR  beam  model 
will  be  generated  next  and  finally,  the  integrated  SOR  beam/shell  model  will  be  generated. 
The  original  material  used  for  all  four  models  is  as  specified  in  Table  1.  Only  the  smearing 
values  will  be  different.  This  section  will  conclude  with  a  comparison  of  the  displacement, 
velocity  and  acceleration  responses  for  the  models  subjected  to  both  a  symmetric  and  off- 
center  shock  wave. 

A.         DISCRETE  STIFFENED  ONR  CYLINDER 

The  three-dimensional  model  is  shown  in  Figure  32.  As  with  the  unstiffened  model, 
a  full  model  is  depicted  even  though  a  half  model  was  analyzed  due  to  symmetry  conditions. 
Both  the  unstiffened  and  stiffened  models  contain  1080  wet  shell  elements,  however  with 
the  stiffeners,  the  total  amount  of  elements  for  the  stiffened  model  is  2064.  The  term  wet 
means  exposed  to  the  shock  wave  pressure.  This  type  of  model  involves  longer 
computational  time  to  complete  the  computation  of  the  response  to  the  shock  wave  due  to 
the  increased  number  of  wet  elements.  It  is  because  of  both  the  computational  time  and 
design  effort  that  a  simpler  model  is  needed. 

The  stiffened  or  discretized  model  was  subjected  to  both  symmetric  and  off-center 
shock  waves.  As  with  the  unstiffened  model,  there  was  no  global  deformation  observed  in 
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Figure  32.  ONR  Full  Cylinder  Model  with  Internal  Stiffeners 

the  x-axis  (see  Figure  32).  The  global  deformation  responses  are  provided  in  Figure  33  at 
discrete  states  of  time.  The  plots  are  again  magnified  at  a  scale  factor  of  50  for  viewer 
enhancement.  The  deformation  analysis  indicates  that  there  is  a  large  local  deformation 
gradient  at  the  locations  of  the  deep  frames  and  endplates.  The  gradient  smoothed  out  as  the 
shock  wave  progresses  through  the  cylinder  and  by  3  msec,  there  is  very  little  pinching  effect 
caused  by  the  deep  frames.  Displacement,  velocity,  and  acceleration  plots  will  be  presented 
later  in  this  section. 
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Figure  33a)    State  time  of  0.5  msec 


Figure  33b)    State  time  of  1.0  msec 


Figure  33c)    State  time  of  1.5  msec 
Figure  33a,b,c.  Deformation  Results  of  Stiffened  ONR  Cylinder  Model 
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Figure  33d)  State  time  of  2.0  msec 


Figure  33e)  State  time  of  2.5  msec 


Figure  33f)  State  time  of  3.0  msec 
Figure  33d,e,f  Deformation  Results  of  Stiffened  ONR  Cylinder  (cont.) 
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B.         ORTHOTROPICALLY  SMEARED  ONR  CYLINDER  MODEL 

The  first  simplification  model  of  the  stiffened  ONR  cylinder  investigated  was  the 
orthotropically  smeared  ONR  cylinder.  This  type  of  simplification  attempts  to  properly 
model  the  cylinder  as  one  unit  with  homogenized  properties.  Once  the  effective  or  smeared 
properties  are  computed,  the  cylinder  can  be  subjected  to  the  underwater  explosion  and  the 
results  analyzed.  The  DYNA3D  code  allows  the  user  to  select  an  orthotropic  elastic 
material.  The  material  selection  requires  9  properties  in  order  to  complete  the  material 
matrix.  The  six  properties  would  represent  the  smeared  properties  of  the  original  ONR 
cylinder.  The  first  three  properties  were  the  Poisson's  ratio  for  all  three  material  axes.  The 
material  axis  corresponds  to  the  axis  shown  in  Figure  1 .  All  three  of  these  values  were 
assumed  to  be  the  same  value  of  the  original  material  or  0.3. 

Next  the  values  for  the  Young's  modulus  were  developed.  The  value  for  smeared 
Young's  modulus  in  the  axial  direction  (y-direction)  is  simply  the  discrete  value  for  the 
original  material  since  ring  stiffeners  do  not  contribute  much  in  the  longitudinal  direction. 
To  determine  the  smeared  elastic  modulus  values  in  the  radial  (x  or  z-direction)  and 
tangential  direction  (x  or  z-direction),  the  assumption  was  made  that  the  radial  stiffness  of 
the  shell  is  significantly  less  important  than  the  tangential  stiffness.  The  tangential  stiffness 
in  other  words  directly  affects  the  bending  stiffness  of  the  cylinder.  This  assumption  is  also 
stated  by  Timoshenko  and  Woinowsky-Kreiger  [Ref.  9].  Therefore  the  smeared  value  in 
the  tangential  direction  would  be  set  equal  to  the  radial  direction  smeared  value  (  Ee.  =  Er ). 
To  calculate  this  smearing  value  for  the  radial  and  tangential  direction,  the  2-dimensional 
static  smearing  technique  developed  in  [Ref.  9]  was  used.     This  smearing  value  would  be 
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used  in  the  three-dimensional  representation  of  the  stiffened  ONR  cylinder.  The  smearing 
value  for  both  directions  was  computed  by  the  following  formula  [Ref.  9]  : 


+  [12(1-  v*)    £^ 

h3  a 


where: 


E'  =  E9=E, 

E=  discrete  or  original  Young's  modulus  for  cylinder 

v=  Poisson's  ratio 

h=  thickness  of  cylinder 

Es  =  Young's  modulus  for  stiffeners 

Is  =  moment  of  inertia  of  the  stiffener 

The  remaining  three  values  are  for  the  elastic  shear  modulii  for  all  three  planes.  The  values 

for  Gxy '  and  G^ '  were  assumed  equal  to  each  other  due  to  the  same  reasoning  for  the  elastic 

modulus.  The  value  for  these  two  planes  was  found  by  using  the  following  formula  [Ref.  9]: 


G'xy=   G'yz  = 


2(l+v) 


where 

E=  elastic  modulus  of  original  cylinder 
v=  Poisson's  ratio  for  original  cylinder 

The  value  for  the  elastic  shear  modulus  G^ '  was  found  using  the  following  expression 

[Ref.9]:: 

G 

GL  -  G  +  -1 
2a 
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where 

G=  original  elastic  shear  modulus  of  cylinder 
Gs  =  shear  modulus  of  stiffeners 
a=  spacing  between  stiffeners 

Besides  ensuring  that  the  proper  stiffness  smeared  values  are  inserted  into  the 

smeared  model,  the  effective  mass  distribution  must  also  be  accounted  for.    The  effective 

mass  distribution  for  the  smeared  model  was  calculated  by  using  the  formula: 


V 


where: 

p'  =  smeared  mass  density  value 

M  =  total  mass  including  cylinder  and  stiffeners  to  be  smeared 

V=  volume  of  cylindrical  shell  without  stiffeners 


Due  to  the  large  thickness  of  the  deep  frame  stiffener,  the  deep  frames  were  discretely 
modeled  in  the  smeared  ONR  cylinder  model.  The  model  is  presented  in  Figure  34  with 
the  axis  as  shown.  Table  3  summarizes  the  inputs  for  the  orthotropic  model.  This  model 
was  then  subjected  to  the  underwater  explosive  testing  for  both  symmetric  and  off-center 
shot  geometries.  The  overall  global  deformation  plots  for  the  symmetric  shot  geometry  is 
provided  in  Figure  35  at  a  magnified  scale  factor  of  50.  There  is  slight  pinching  around  deep 
frames  at  the  earlier  state  times  but  the  deformation  gradient  around  each  of  the  deep  frames 
is  less  than  the  original  discrete  model.  As  will  be  seen  in  futures  ection  in  this  chapter,  the 
overall  dynamic  reponses  of  this  model  are  in  close  agreement  with  the  discrete  model. 
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Figure  34.    Orthotropically  Smeared  Model  of  ONR  Stiffened  Cylinder 


v(xz,yz,xy) 

0.3 

0.3 

0.3 

E  (psi) 

(Exx>  Eyy>  Ezz) 

17.813  x  108 

same 

2.9  xlO6 

G(psi) 

(xz,  yz,  xy) 

12.145  x  106 

11.154  x  106 

11.154  xlO6 

Table  3.  Material  Properties  for  Orthotropic  Smeared  ONR  Cylinder 
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The  use  of  the  orthotropic  smeared  technique  produces  results  similar  to  the  study  on 
the  flat  plate  where  local  smearing  was  used  around  the  clamped  edges.  In  both  cases,  local 
smearing  reduced  the  deformation  gradient  and  produced  effective  global  deformation 
results.  One  of  the  disadvantages  of  using  the  orthotropic  smeared  model  is  more  lengthy 
computation  times  due  to  explicit  time  integration  scheme.  Although  the  model  was 
reduced  from  2064  nodes  to  1 124  nodes,  the  amount  of  wet  elements  remains  the  same. 
Since  smearing  also  involves  increasing  the  Young's  modulus,  the  overall  critical  time  step 
is  decreased  from  the  original  isotriopic  model's  time  step.  These  factors  cause  a  very 
lengthy  computational  time. 
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Figure  35a)    State  time  of  0.5  msec 


Figure  35b)    State  time  of  1.0  msec 


Figure  35c)    State  time  of  1.5  msec 
Figure  35a,b,c.  Deformation  Results  of  Orthotropically  Smeared  ONR  Cylinder  Model 
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Figure  35d)  State  time  of  2.0  msec 


Figure  35e)  State  time  of  2.5  msec 


Figure  35f)  State  time  of  3.0  msec 
Figure  35d,e,f.  Deformation  Results  of  Orthotropically  Smeared  ONR  Cylinder  (conL) 
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C.         SOR  BEAM  MODEL 

The  SOR  beam  model  for  the  stiffened  ONR  cylinder  is  developed  in  the  same 
procedure  as  the  unstiffened  SOR  beam  model.  However,  the  effective  stiffness  and  mass 
distribution  of  the  stiffened  model  must  be  incorporated  into  the  SOR  beam  model.  This  will 
ensure  that  the  beam  model  will  have  the  equivalent  fluid  structure  interaction  as  compared 
with  the  stiffened  ONR  cylinder.  The  same  harmonic  function  will  be  used  in  computing  the 
fluid  mass  matrix  for  the  SOR  beam  model.  Both  symmetric  and  off-center  shot  geometry 
will  be  used  for  analyzing  the  model. 

There  are  45  SOR  beam  elements  used  to  model  the  cylinder  along  its  longitudinal 
direction.  In  order  to  obtain  the  equivalent  stiffness,  the  cross-sectional  areas  and  moments 
of  inertia  were  adjusted.  The  cross-section  area  of  the  stiffened  SOR  beam  model  is: 


A  -  ±tf  -  df) 
4 


where  d0  and  d^  are  the  outside  and  inside  diameters  respectively.  This  value  is  the  same  for 
both  stiffened  and  unstiffened  models.  However,  the  second  moment  of  inertia  and  polar 
moment  of  inertia  of  the  stiffened  cylinder  differed  from  the  unstiffened  cylinder  because 
the  stiffeners  were  considered.    The  values  were  computed  as  follows: 


I  '    [~7A    °     "      i  "cylinder   +    L~("o    "   "/  'ideepframe   +    L  ~  ("o    "    «J  n  shallow  frame 
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J  =  2*  / 

The  mass  distribution  value  was  the  same  as  the  orthotropically  smeared  mass  value. 
Figure  36  shows  the  SOR  beam  model  of  the  stiffened  ONR  cylinder.  The  model  was 
subjected  to  shock  wave  and  the  deformation  plots  were  made  and  shown  in  Figure  37.  The 
deformation  has  a  similar  shape  to  the  discrete  model,  however  closer  examination  revealed 
that  the  deformation  was  in  fact  less  than  the  discrete  model.  This  was  due  to  not  taking  into 
account  the  local  effects  of  the  deep  frames  and  endplates.  Since  the  mass  distribution  and 
cross-sectional  properties  of  the  model  only  accounted  for  the  shallow  frames  ,  the  model 
doe  not  properly  show  the  local  mass  density  and  equivalent  stiffness  of  the  endplates  and 
deep  frames. 
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Figure  36.  SOR  Beam  Model  of  Stiffened  ONR  Cylinder 

63 


y   < 


i    I    '    i    i    '    f    I    I    I    '    I    I    M    I    '    [    L 


Figure  37a)    State  time  of  0.5  msec 
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Figure  37b)    State  time  of  1.0  msec 
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Figure  37c)    State  time  of  1.5  msec 


Figure  37a,b,c.  Deformation  Results  of  SOR  Beam  Model  of  ONR  Cylinder  Model 
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Figure  37d)  State  time  of  2.0  msec 
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Figure  37e)  State  time  of  2.5  msec 
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Figure  37f)  State  time  of  3.0  msec 
Figure  37  d,e,f.  Deformation  Results  for  SOR  Beam  Model  of  ONR  Cylinder  (cont) 
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To  correct  the  mass  distribution  and  equivalent  stiffness  of  the  SOR  beam  model, 
the  mass  density  of  the  deep  frames  and  endplates  were  inputted  at  their  respective  element 
location  in  the  DYNA3D  program.  The  cross  sectional  properties  were  also  adjusted 
locally  for  the  deep  frames  and  endplates.  Additionally,  the  element  length  for  the  deep 
frame  section  was  reduced  to  two  inches  which  corresponds  with  the  original  cylinder 
dimensions.  This  allowed  for  a  more  accurate  local  smearing  of  both  mass  density  and  cross- 
sectional  properties  for  the  deep  frames.  The  adjusted  SOR  Beam  model  now  contains  49 
elements  and  is  shown  in  Figure  38.  The  global  deformation  plots  are  presented  in  Figure 
39.  The  deformation  plots  show  an  improvement  from  the  globally  mass  smeared  45  SOR 
beam  element  model.  Figures  40-42  show  the  axial  dynamic  responses  at  the  left  endplate 
which  support  using  local  smearing  of  the  mass  density  for  deep  frames  and  endplates. 
Detailed  comparisons  between  discrete  and  SOR  beam  model  will  be  presented  later  in  this 
report.    The  computational  time  using  the  SOR  beam  model  is  approximately  1/20  of  the 

discrete  model's  computational  time. 
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Figure  38.  SOR  Beam  Model  With  Local  Mass  Density  Smearing 
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Figure  39a)  State  time  of  0.5  msec 
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Figure  39b)  State  time  of  1.0  msec 
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Figure  39c)  State  time  of  1.56  msec 


Figure  39a,b,c.  Deformation  Results  For  SOR  Beam  Model  With  Local  Mass  Density 
Smearing  of  Endplates  and  Deep  Frames 
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Figure  39d)  State  time  of  2.0  msec 


Figure  39e)  State  time  of  2.5  msec 
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Figure  39f)  State  time  of  3.0  msec 


Figure  39d,e,f.  Deformation  Results  for  SOR  Beam  Model  With  Local  Mass  Density 
Smearing  of  Deep  Frames  and  Endplates  (cont.) 
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Figure  40.    Axial  Displacement  Response  Comparing  Local  to  Global  Mass  Density 
Smearing  of  SOR  Beam  Model 
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Figure  41.    Axial  Velocity  Response  Comparing  Local  to  Global  Mass  Density 
Smearing  of  SOR  Beam  Model 
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Figure  42.    Axial  Acceleration  Response  Comparing  Local  to  Global  Mass  Density 
Smearing  of  SOR  Beam  Model 
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D.        INTEGRATED  SOR  BEAM/SHELL  MODEL 

The  integrated  SOR  beam/shell  model  incorporates  all  simplification  techniques 
discussed  in  this  paper.  The  interface  will  be  modeled  using  the  rigid  body  material  option 
and  local  smearing  will  be  done  for  the  beam  elements  to  take  into  account  the  mass 
distribution  of  the  deep  frames  and  endplates.  Since  it  was  desired  to  examine  the  dynamic 
responses  at  one  of  the  deep  frames,  a  beam/shell  model  was  produced  with  fewer  shell 
elements  as  shown  in  Figure  43.  The  model  contains  a  total  of  450  shell  elements  and  29 
SOR  beam  elements.  The  shell  section  has  the  shallow  frames  and  deep  frames  discretely 
modeled  to  produce  the  fastest  computational  time.  The  beam  section  has  the  cross-section 
and  second  moment  of  inertia  values  adjusted  to  account  for  the  shallow  and  deep  frame 
stiffeners.  This  model  will  be  tested  with  a  charge  located  off-center  as  shown  in  Figure 
44.  The  symmetric  test  on  beam/shell  model  would  have  only  produced  a  larger  model  and 
hence  it  was  not  tested. 


x 
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Figure  43.  SOR  Beam/Shell  Model  of  Stiffened  ONR  Cylinder 
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25"  from  center 
of  cylinder 


HBX-1  charge 


Figure  44.  Offcenter  Shot  Geometry  for  Integrated  Beam/shell  Model 

The  global  deformation  plots  are  provided  in  Figure  45.  The  deformation  plot  shows 
local  pinching  at  deep  frames  and  endplates  which  is  consistent  with  the  discrete  model.  The 
rigid  body  model  of  the  dry  endplate  produces  larger  local  deformation  at  the  endplate.  The 
overall  global  deformation  appears  to  be  in  good  agreement  with  the  other  three  models. 

The  computational  time  for  the  model  was  about  1/3  of  the  discrete  model  which 
clearly  supports  using  this  type  of  model.  It  is  highly  recommended  to  use  as  many  SOR 
Beam  elements  as  possible  in  the  model  design  to  reduce  the  computational  time.  Detailed 
analysis  at  various  locations  on  the  beam/shell  model  will  be  provided  later  in  this  chapter. 
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Figure  45a)   State  time  of  0.5  msec 


Figure  45b)    State  time  of  1.0  msec 


Figure  45c)    State  time  of  1.5  msec 


Figure  45a,b,c.  Deformation  Results  of  Integrated  SOR  Beam/Shell  Model  of  ONR 
Cylinder  Model 
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Figure  45d)  State  time  of  2.0  msec 


Figure  45e)  State  time  of  2.5  msec 


Figure  45f)  State  time  of  3.0  msec 
Figure  45d,e,f  Deformation  Results  for  Integrated  SOR  Beam/Shell  Model  (cont.) 
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E.         COMPARISON  OF  DISCRETE  TO  SIMPLIFIED  MODELS 

This  section  will  compare  the  dynamic  responses  between  the  fully  discretized 
stiffened  ONR  cylinder  to  the  simplified  models.  The  analysis  will  compare  both  symmetric 
and  off-center  geometry  for  the  location  of  the  charge.  The  displacement,  velocity,  and 
acceleration  responses  will  be  displayed  and  discussed. 

1.         Comparison  Using  Symmetric  Explosive  Testing 

The  geometry  for  symmetric  testing  is  the  same  as  Figure  2.  The  aim  of  this  testing 
was  to  simply  compare  the  responses  between  the  orthotropically  smeared,  beam,  and 
discrete  models  of  the  stiffened  ONR  cylinder.  The  beam/shell  model  was  not  included  since 
more  detailed  data  was  taken  using  off-center  shot  geometry.  The  results  will  aim  to 
compare  the  responses  at  a  deep  frame  location,  endplates  and  at  the  shot  location. 

Figure  46  compares  the  overall  global  longitudinal  deformation  plots  for  all  four 
models  at  state  time  0.5  msec  and  at  magnified  scale  factor  of  50.  All  plots  are  taken  at  the 
same  coordinate  axis  position.  There  is  good  agreement  except  for  areas  of  local 
deformation  at  the  endplates  and  deep  frames.  The  deformation  plots  show  that  both  the 
smeared  stiff ener  technique  and  SOR  beam  element  modeling  are  effective  in  producing 
similar  global  responses.  However,  it  is  recommended  that  use  of  large  numbers  of 
orthotropic  shell  elements  should  be  avoided  since  longer  computational  times  occur  when 
using  these  types  of  elements  due  to  the  decreased  critical  time  step. 
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Figure  46a)    ONR  Stiffened  Cylinder  at  0.5  msec 


Figure  46b)    Orthotropically  Smeared  ONR  Cylinder  at  0.5  msec 
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Figure  46c)    SOR  Beam  Model  of  ONR  Cylinder  at  0.5  msec 


Figure  46a,b,c  Deformation  Results  for  Discrete  vs.  Simplified  Models  for  ONR 
Stiffened  Cylinder  with  Symmetric  Underwater  Explosion  Testing 
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Next,  the  axial  displacement,  velocity,  and  acceleration  plots  were  made  at  the  left 
and  right  endplate  areas  to  check  global  motion  at  the  ends.  Figures  47-49  shows  the 
displacement,  acceleration  of  all  four  models  at  the  left  endplate.  There  is  a  shift  in  the 
frequency  between  the  SOR  beam  and  discrete  model,  however  the  beam's  peak  magnitude 
is  very  close  to  the  shell's  first  oscillation  peak  magnitude.  There  is  significantly  less 
vibration  and  acceleration  experienced  in  the  beam  model  as  compared  to  the  shell  model 
which  is  expected.  The  orthotropic  model  compares  favorably  to  the  discrete  model  in  shape 
but  is  slightly  off  in  magnitude. 
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Figure  47.  Axial  Displacement  Comparison  at  Left  Endplate 
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Figure  48.  Axial  Velocity  Comparison  At  Left  Endplate 
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Figure  49.  Axial  Acceleration  Comparison  at  Left  Endplate 
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Figure  50-52  shows  the  whipping  responses  for  all  three  models  at  the  shock  point 
(nearest  structure  point  from  the  charge).  The  orthotropic  shell  model  and  discrete  model 
are  very  similar  in  magnitude  and  shape  and  have  an  average  error  of  less  than  15%.  There 
is  less  velocity  and  acceleration  in  the  SOR  Beam  model  because  the  local  smearing 
represents  the  average  stiffness  which  produce  fewer  oscillations  than  the  original  model. 
The  cylindrical  models  have  radial  deformation  which  produce  more  velocity  and 
acceleration  responses  as  the  shock  wave  progresses  through  the  shell. 
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Figure  50.  Whipping  Displacement  Comparison  at  Shock  Point 
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Z-Direction  Velocity  for  All  Three  ONR  Models 
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Figure  51.  Whipping  Velocity  Comparison  at  Shock  Point  of  All  Three  Models 
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Figure  52.  Whipping  Acceleration  Comparison  at  Shock  Point  For  All  Three  Models 
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And  finally,  the  dynamic  responses  are  analyzed  at  a  deep  frame  on  the  three  models.  This 
location  is  75  inches  from  the  center  of  the  cylinder  (see  Figure  1).  Only  the  whipping 
motion  responses  will  be  presented,  however  there  was  good  agreement  in  the  axial  direction 
at  the  deep  frame.  There  was  reduced  velocity  and  acceleration  in  axial  direction  which  was 
consistent  with  the  endplate  responses.  Figures  53-55  show  the  whipping  responses  at  the 
deep  frame.  The  results  show  excellent  agreement  between  the  discrete  and  SOR  beam 
model  with  only  slight  differences  in  velocity  and  acceleartion  for  reasons  previously  noted. 
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Figure  53.  Whipping  Displacement  Comparison  at  Deep  Frame 
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Figure  54.  Whipping  Velocity  Comparison  at  Deep  Frame  Location 
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Figure  55.  Whipping  Acceleration  Comparison  At  Deep  Frame 
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Overall,  the  dynamic  responses  can  support  using  either  SOR  beam  elements,  or 
orthotropic  shell  elements.  There  are  tradeoffs  to  be  considered  for  choosing  which 
technique  to  use.  The  orthotropic  shell  element  produces  more  accurate  responses  than  the 
SOR  beam  element,  however  the  SOR  beam  element's  computational  time  is  about  20  times 
faster  than  an  orthotropic  shell  elements.  When  only  49  elements  are  used  for  SOR  beam 
model  vice  1080  for  the  orthotropic  shell  model,  there  is  a  big  reduction  in  computational 
time,  but  with  the  reduced  accuracy.  So  both  techniques  should  be  considered  for 
analyzing,  then  if  time  is  the  more  important  factor  in  the  solution,  use  of  SOR  beam 
elements  are  recommended.  Table  4  provides  a  comparison  based  on  size  of  model  and 
overall  computational  time  for  all  three  models  subjected  to  symmetrci  shock  testing. 


Stiffened  ONR 
Cylinder 

SOR  Beam  Model 

Orthotropic 
Smeared  Model 

Number  of 
Nodes 

2064 

51 

1224 

Number  of  Wet 
Elements 

1080 

49 

1080 

Computational 
Time  (sec) 

2070 

95 

12200 

Table  4.  Comparison  Data  For  all  Three  Models  of  Stiffened  ONR  Cylinder 
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2.  Comparison  Using  Off-center  Explosive  Testing 

This  section  will  compare  the  dynamic  responses  between  the  fully  discretized 
stiffened  ONR  cylinder  to  the  SOR  Beam  and  integrated  beam/shell  model  for  an  off-center 
shot  geometry.  The  ONR  stiffened  model  was  subjected  to  the  underwater  explosion  at  a 
location  25  inches  to  the  right  of  the  center  of  the  cylinder  (See  Figures  1  &44).  The  SOR 
beam  model  was  the  same  model  used  in  the  symmetric  analysis  but  subjected  to  the  off- 
center  shock  wave.  Finally,  the  integrated  beam/shell  model  was  subjected  to  the  same 
testing.  The  integrated  beam/shell  model  used  the  rigid  body  material  option  for  the  endplate 
and  locally  smeared  the  mass  distribution  and  cross-sectional  properties  of  the  deep  frames 
and  endplates.  The  orthotropic  shell  model  was  not  considered  due  to  its  lengthy 
computational  time. 

Figure  56  is  a  global  deformation  plot  of  all  three  models  at  0.5  msec.  The  models 
appear  to  be  in  good  agreement  except  for  some  local  pinching  at  the  dry  endplate  on  the 
integrated  beam/shell  model.  There  is  however  some  difference  in  shape  between  the  SOR 
beam  model  and  discrete  model  which  is  mainly  due  to  applying  the  average  stiffness  and 
mass  distribution  at  the  deep  frames  and  endplates.  This  approximation  produces  similar 
results  but  will  not  exactly  match  the  discrete  model. 

Since  the  main  objective  in  this  study  was  to  simplify  the  original  model,  it  is 
important  to  select  an  area  on  the  shell  section  of  the  integrated  beam/shell  model  and 
compare  the  responses  between  discrete  and  the  simplified  model.  If  the  responses  are  in 
good  agreement  at  this  selected  area,  then  further  internal  structures  could  be  placed  inside 
the  shell  section  of  the  beam/shell  model.  Thus,  only  a  small  portion  of  the  original  model 
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Figure  56a)    ONR  Stiffened  Cylinder  at  0.5  msec 
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Figure  56b)    SOR  Beam  Model  of  ONR  Cylinder  at  0.5  msec 


Figure  56c)    Integrated  SOR  Beam/Shell  Model  of  ONR  Cylinder  at  0.5  msec 


Figure  56a,b,c.  Global  Deformation  Plots  -  Off-center  Explosion  Testing  for  Stiffened 
ONR  Cylinder 
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would  need  to  be  modeled  discretely,  and  the  rest  of  the  model  could  be  modeled  with  SOR 
beam  elements.  The  selected  area  for  comparison  purposes  was  the  deep  frames  which  is 
located  75  inches  right  of  the  center  of  the  original  cylinder  and  is  discretely  modeled  inside 
the  shell  section  of  the  beam/shell  model.  Both  axial  and  whipping  responses  were 
analyzed.  Figures  57-59  show  the  axial  responses  and  Figure  60-62  show  the  whipping 
responses  at  the  deep  frame.  Since  the  integrated  beam/shell  model  used  SOR  beam 
elements  for  60%  of  the  longitudinal  length  of  the  original  cylinder,  slightly  higher 
inaccuracies  were  experienced.  The  integrated  beam/shell  model  responses  are  in  closer 
agreement  to  the  discrete  model.  The  whipping  responses  are  very  similar  in  shape  and 
only  differ  in  magnitude  due  to  the  use  of  SOR  beam  elements.  Table  5  is  provided  for 
comparing  the  size  and  computational  time  for  all  three  models. 


Stiffened  ONR 
Cylinder 

SOR  Beam  Model 

Integrated 
Beam/shell  Model 

Number  of 
Nodes 

2064 

51 

942 

Number  of  Wet 
Elements 

1080 

49 

479                        i 

Computational 
Time  (sec) 

2070 

95 

719 

Table  5.  Comparison  Data  For  all  Three  Models  of  Stiffened  ONR  Cylinder 


86 


0.14 


Axial-Displacement  for  All  Three  ONR  Models 


1  1.5  2 

Time,  ( milli  seconds) 


Figure  57.  Axial  Displacement  Comparison  At  Deep  Frame  Location 


87 


200 


150 


100 


50 


Axial-Velocity  for  All  Three  ONR  Models 


-50 


-100 


SOR  Beam 


ONR  Orig  Model 
Integrated  Beam/Shell  Model 


0.5 


1  1.5  2 

Time,  ( milli  seconds) 


2.5 


Figure  58.  Axial  Velocity  Comparison  At  Deep  Frame  Location 
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Figure  59.  Axial  Acceleration  Comparison  At  Deep  Frame  Location 
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Z-Direction  Displacement  for  all  Three  ONR  Models 
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Figure  60.  Whipping  Motion  Comparison  At  Deep  Frame  Location 
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Figure  61.  Whipping  Motion  Comparison  At  Deep  Frame  Location 


Z-Direction  Acceleration  for  All  Three  ONR  Models 
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Figure  62.  Whipping  Motion  Comparison  At  Deep  Frame  Location 
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Overall,  the  responses  favor  using  an  integrated  beam/shell  model  if  additional 
internal  structures  require  modeling.  This  model  provides  both  quicker  computational  time 
and  required  accuracy  of  dynamic  responses. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

Simplified  finite  element  modeling  of  a  stiffened  structure  is  an  effective  means  by 
which  both  the  design  effort  and  the  computational  time  can  be  reduced.  When  the  discrete 
and  simplified  models  are  subjected  to  underwater  explosions  and  closely  examined,  the 
results  show  that  either  smearing  orthotropically  or  using  the  SOR  beam  element  produces 
similar  global  results.  When  examining  the  two  models  at  very  discrete  points,  the  results 
are  similar  but  differ.  The  orthotropically  smeared  cylindrical  model  produces  more  accurate 
comparisons,  but  has  a  very  lengthy  computational  time.  There  is  more  inaccuracy  using 
SOR  beam  element  modeling,  however  the  computational  time  and  design  effort  is 
significantly  reduced. 

Incorporating  the  advantages  of  both  techniques  into  an  integrated  beam/shell  model 
provides  the  best  means  of  simplifying  the  discrete  model.  The  results  are  more  favorable 
and  the  computational  time  is  still  significantly  reduced.  The  interface  between  the  beam 
and  shell  section  can  be  modeled  as  either  a  dry  endplate  with  smeared  properties  or  by  using 
the  rigid  body  material  option  within  DYNA3D  code.  It  is  recommended  that  the  rigid  body 
material  option  be  used  since  the  critical  time  step  is  not  changed  when  this  material  is 
selected  for  modeling  the  interface.  The  use  of  the  beam/shell  model  allows  further  internal 
structures  to  be  modeled  without  dramatically  increasing  both  design  effort  and 
computational  time  for  numerical  solution.  If  a  quick  but  less  accurate  solution  is  needed 
for  a  particular  problem,  then  the  SOR  beam  model  is  recommended. 
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Since  no  physical  experimental  data  was  used  in  this  study,  it  is  highly  recommended 
that  a  scaled  down  experimental  model  be  built  and  tested  against  the  numerical  data 
collected  in  this  study.  The  present  numerical  analysis  involved  using  an  explicit  time 
integration  scheme  adopted  in  DYNA3D  code.  An  additional  comparison  of  the  results 
could  be  accomplished  by  using  an  unconditionally  stable  implicit  time  integration 
technique  with  another  finite  element  program. 
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APPENDIX  A 

INPUT  FILES  FOR  SOR  BEAM  MODEL 
The  following  input  files  consist  of  ingrid  which  creates  the  geometric  model 
for  DYNA3D,  flumas,  augmat,  and  timint.    The  flumas,  augmat  and  timint  files 
provide  input  for  the  USA  program.  These  files  are  for  SOR  Beam  Model  Only. 
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35 

1 

36 

1 

37 

1 

38 

1 

39 

1 

40 

2 

41 

1 

42 

1 

43 

1 

44 

1 

45 

1 

46 

1 

47 

1 

48 

1 

49 

3 

25 
26 

27 
28 
29 

30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 


26 

51 

27 

51 

28 

51 

29 

51 

30 

51 

31 

51 

32 

51 

33 

51 

34 

51 

35 

51 

36 

51 

37 

51 

38 

51 

39 

51 

40 

51 

41 

51 

42 

51 

43 

51 

44 

51 

45 

51 

46 

51 

47 

51 

48 

51 

49 

51 

50 

51 

LOAD  CURVE  DEFINITIONS 


1     2 
.OOOE+OO  2.000E-06 
1.000E+00  2.000E-06 


FLOMAS  DATA  FOR  MODEL  BEAMONR 

beam. flu  beam . geo  strnam.dat  beam.daa 


T  T  T  T 

T  F  F  F 

F  F  F  T 

F  F  F  F 

F  T  F  F 

F  F  F  F 

F  F  F  F 

DYNA 

0  51  0  0 

10   0 

0   112 

0 

51 

0 

9.356E-05 

57600 

.0 

5 

1    50   16 

.8   16 

.  8 

3 

1     10 

16 

.8 

1 

1    50   16 

.8   16 

.6 

49 

50   50   16 

.8    0 

1 

$   FLUNAM   GEONAM   GRDNAM   DAANAM 

$ 

PRTGMT 

PRTTRN 

PRTAMF 

CALCAM 

$ 

EIGMAF 

TWODIM 

HAFMOD 

QUAMOD 

$ 

PCHCDS 

NASTAM 

STOMAS 

STOINV 

$ 

FRWTFL 

FRWTGE 

FRWTGR 

FRESUR 

$ 

RENUMB 

STOGMT 

ROTGEO 

ROTQUA 

$ 

PRTCOE 

STRMAS 

SPHERE 

ROTSYM 

$ 

OCTMOD 

CAVFLU 

FRWTFV 

INTCAV 

$ 

MAINKY 

$ 

NSTRC 

NSTRF 

NGEN 

NGENF 

$ 

NBRA 

NCYL 

NCAV 

$ 

NHAS 

NHAF 

NHAI 

NFUN 

$ 

NSEG 

$ 

NCIR 

$ 

RHO 

CEE 

S 

NVEC 

s 

Nl 

N2 

Rl     R2 

NSET 
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AUGMAT    DATA    FOR    MODEL . BEAMONR 
strnam.dat    beam.daa   beam . geo    beam   pre 
F    F    F   F 


T    F 

F    F 
F    F 


F 

F 

F 

F 

DYNA 

0.0 

51    1S3 


3    3 


S       STRNAM       FLUNAM      GEONAM       PRENAM 
$      FRWTGE      FRWTST      FRWTFL      PLNWAV 
FLUSKY      DAAFRM      SYMCON      DOFTAB 
PRTGMT      PRTTRN      PRTSTF      PRTAUG 
MODTRN      STRLCL       INTWAT      CFADYN 

MAINKY 
NTYPDA 


r=d3dump 

TIME    INTEGRATION 

DATA 

FOR    MODEL    ! 

BEAMONR 

beam. pre    beam.pos 

$ 

PRENAM 

POSNAM 

beam. rst 

$ 

RESNAM 

F    F    F 

$ 

BUBPUL 

REFSEC 

F LUMEN 

T    F    F    F 

$ 

EXPWAV 

SPLINE 

VARLIN 

PACKET 

F   T    F    F 

$ 

HYPERB 

EXP LOS 

DOUBDC 

VELIMP 

1 

$ 

NCHARG 

0.0 

$ 

HYDPRE 

0.0    -25 .3125    326  . 

4 

$ 

XC         YC 

ZC 

0.0     -25.3125       16. 

813 

$ 

SX         SY 

SZ 

101 

$ 

JPHIST 

1.0 

$ 

PNORM 

14.0E-06 

$ 

DTHIST 

1 

$ 

CHGTYP 

60.0       25.8       10.0 

$ 

WEIGHT 

SLANT 

CHGDEP 

1    0 

$ 

NTINT 

NCHGAL 

0.0       2.00E-06 

$ 

STRTIM 

DELTIM 

5000       5000 

$ 

NSAVER 

NRESET 

0    0    0    0 

$ 

LOCBEG 

LOGRES 

LOCWRT 

NSTART 

F    F 

$ 

FORWRT 

STBDA2 

F 

$ 

DISPLA 
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APPENDIX  B 

INPUT  FILES  FOR  SOR  BEAM/SHELL  MODEL 
The  following  input  files  are  for  the  integrated  SOR  Beam/shell  model. 


bsmodel .half 

9  942  0  32  860 


88    large 


ooooo 

0                         0  e20.9  450         33 

2             0       450              0000             0  0000              0  000 

3.000E-03     1.000E+03     1.000E-0S  .OOOE+00  .OOOE+00              0  .000E+00 

0000020010  0000 
0            0        .OOOE+00 

0  0  0 

1  37.3560E-04            0           0  .OOOOE+00  0     .OOOOE+00    .00O0E+OO  2            0           0 
material    type    #    3        (elastic-plastic) 

2.900E+07        .0OOE+00        .OOOE+00  .000E  +  00  .OOOE  +  00  OOOE  +  00        -OOOE  +  00        .  000E  +  00 

3.000E-01        .OOOE+00       .0O0E+00  .000E+00  .000E+00  OOOE+00        .OOOE+00       .OOOE+00 

2.000E+10       .0OOE+00       .OOOE+00  .000E+00  .000E+00  OOOE+00        .OOOE+00        .000E+00 

5.020E+04        .OOOE+00       .OOOE+00  .000E+00  .000E+00  .OOOE+00        .OOOE+00        .000E+00 

.OOOE+00        .00OE+00        .OOOE+00  .OOOE+00  .OOOE+00  .OOOE+00        .OOOE+00          OOOE+00 

OOOE+00        .000E+00        .OOOE+00  .000E+00  . OOOE+00  .OOOE+00        .OOOE+00        .000E+0C 
section    properties 

5 

1.880E-01    1.88OE-01    1.880E-01  1.880E-01  . OOOE+00 

2  17.3560E-04            0           0  .0000E  +  00  0     .OOOOE+00      OOOOE+OO  10            2 
material    type    8    1       (elastic) 

2.900E+07        .OOOE+00       .000E+00  .000E+00  .OOOE+00  .OOOE+00        .OOOE+00        .000E+00 

3.000E-01        .00OE+00       .OOOE+00  .OOOE+00  .000E+00  OOOE+00        .OOOE+00       .OOOE+00 

000E+00        .OOOE+00       .OOOE+00  .000E+00  OOOE+00  OOOE+00        .OOOE+00        .OOOE+00 

.OOOE+00       .OOOE+00       .OOOE+00  .000E+00  .000E+00  OOOE+00        .OOOE+00        .OOOE+00 

.OOOE+00       .OOOE+00-     .OOOE+00  .OOOE+00  .OOOE+00  OOOE+00        .OOOE+00        .OOOE+00 

.OOOE+00       .OOOE+00       .OOOE+00  .OOOE+00  .OOOE+00  OOOE+00        .OOOE+00        .OOOE+00 

section    properties 

0 

1.970E+01    5.402E+03    5.402E+03  1.081E+04  1.970E+01  OOOE+00 

3  209.5710E-04            0           0  .0000E+00  0     .0000E+00    .OOOOE+OO  2            0            C 
material    type    #    20       (rigid  body) 

2.900E+07       .OOOE+OO       .OOOE+00  .000E+00  .000E+00  .OOOE+00        .OOOE+00        .OOOE+00 

3.000E-01        .OOOE+OO       .OOOE+00  .OOOE+OO  .OOOE+00  OOOE+00        .OOOE+OO        .OOOE+00 

.OOOE+00       .OOOE+00       .OOOE+OO  .OOOE+OO  .OOOE+00  OOOE+00        .OOOE+00        .OOOE+00 

.OOOE+00       .OOOE+00       .OOOE+OO  .OOOE+OO  .OOOE+OO  .OOOE+00        .OOOE+00        .OOOE+OO 

.OOOE+00       .OOOE+00       .OOOE+00  .OOOE+OO  .OOOE+00  OOOE+00        .OOOE+00        .OOOE+00 

.OOOE+00        .OOOE+00        .OOOE+00  .OOOE+OO  .OOOE+OO  OOOE+00        .OOOE+OO        .OOOE+OO 

section   properties 
5 

1.880E-01    1.880E-01    1.880E-01  1.880E-01  .OOOE+00 

4  37.356OE-04            0           0  .OOOOE+OO  0     .OOOOE  +  00      OOOOE  +  OO  2            0            0 
material    type    #    3       (elastic-plastic) 

2.900E+07       .OOOE+00       .OOOE+00  .OOOE+00  .OOOE+00  .OOOE+OO        .OOOE+00        .OOOE+00 

3.000E-01       .OOOE+00       .OOOE+00  .OOOE+00  .OOOE+OO  .OOOE+00        .OOOE+OO        .OOOE+00 

2.000E+05       .OOOE+00       .OOOE+OO  .OOOE+OO  .OOOE+00  .OOOE+00        .OOOE+00        .OOOE+00 

5.100E+04        .OOOE+00       .OOOE+00  .OOOE+OO  .OOOE+OO  .OOOE+00        .OOOE+OO        .OOOE+OO 

.OOOE+00       .OOOE+00       .OOOE+00  .OOOE+00  .OOOE+00  .OOOE+00        .OOOE+00        .OOOE+00 

.OOOE+OO       .OOOE+00       .OOOE+00  .OOOE+OO  .OOOE+OO  .OOOE+00        .OOOE+OO        .OOOE+OO 
section   properties 

5 

2.500E+00    2.500E+00    2.S00E+00  2.S00E+00  .OOOE+00 

5  37.3560E-04            0           0  .OOOOE  +  OO  0     .OOOOE  +  00    .OOOOE+OO  2            0            0 
material    type    8    3        (elastic-plastic) 

2.900E+07       .OOOE+00       .OOOE+00  .OOOE+OO  .OOOE+00  OOOE+00        .OOOE+00        .OOOE+OO 

3.000E-01       .OOOE+00       .OOOE+00  .OOOE+OO  .OOOE+OO  OOOE+OO        .OOOE+OO        .OOOE+00 

2.000E+05       .OOOE+00       .OOOE+00  .OOOE+00  .OOOE+00  .OOOE+00        .OOOE+00        .OOOE+OO 

5.100E+04        .OOOE+00        .OOOE+00  .OOOE+00  .OOOE+00  .OOOE+OO        .OOOE+00        .OOOE+00 

.OOOE+00        .OOOE+OO        .OOOE+OO  .OOOE+OO  .OOOE+OO  OOOE+00        .OOOE+OO        .OOOE+00 
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.OOOE 
section 

3 .250E+ 
6 
material 
2 . 900E+ 


00  .0 
prope 

00  3  .2 
37.35 
type 

07   .0 


3 .000E 
2 .000E+ 
5 .100E+ 
.000E+ 
.000E+ 
section 

2 .000E+ 

7 

material 

2 .900E+ 

3 .000E- 

2 .000E+ 

5. 100E  + 

.000E 

.000E+ 

section 

2 .500E- 

8 

material 

2 . 900E  + 

3 .000E- 
.000E+ 
.000E+ 
.O0OE+ 
.000E+ 

section 

3.970E+ 

9 

material 

2.900E+ 

3.000E- 
.000E+ 
.000E 
.000E+ 
.000E+ 

section 

8.880E+ 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 


prope 

00  2  .0 
37.35 
type 

07   .0 


prope 

01  2  .5 
18.26 
type 

07   .0 


prope 

02  3  .6 
11.04 
type 

07   .0 


prope 

02  6  .2 
1. 
0. 
0. 
0. 
0. 
0. 

1 . 

0. 
0. 

0. 
0. 

0. 


OOE+OO 
rties 

5 

50E+00 
60E-04 
#3  (e 
00E+00 
00E+00 
00E+00 
00E+00 
00E+00 
00E+OO 
rties 

5 
00E+00 
60E-O4 
#3  (e 
00E+00 
00E+00 
00E+00 
00E+00 
00E+00 
OOE+OO 
rties 

5 
00E-01 
92E-04 

#  1  (e 
OOE+OO 
OOE+OO 
OOE+OO 
OOE+OO 
OOE+OO 
OOE+OO 
rties 

0 
14E+03 
70E-03 

#  1  (e 
OOE+OO 
OOE+OO 
OOE+OO 
OOE+OO 
OOE+OO 
OOE+OO 
rties 

0 
75E+04 
.00000 
.26301 
.51955 
.76329 
.98824 
.11888 
.00000 
.26301 
.51955 
.76329 
.98824 
. 11888 


OOOE  +  OO   .OOOE  +  00   .O00E+00   .OOOE+OO   .  OOOE+OO   . OOOE+OO 


3 .250E+00  3.250E+ 

0     0  .0000E+ 

lastic-plastic) 

.OOOE+00  .000E+ 

.OOOE+OO  .000E+ 

.OOOE+00  .000E+ 

.OOOE+OO  .000E+ 

.OOOE+OO  .000E+ 

.OOOE+00  .000E+ 


2. OOOE+00  2.000E+ 

0     0  .0000E+ 

lastic-plastic) 
.OOOE+OO  .000E+ 
.OOOE+00  .000E+ 
.OOOE+OO  .000E+ 
.OOOE+00  .000E+ 
.OOOE+OO  .00OE+ 
.OOOE+00   .0OOE+ 


2 . 500E-01 
0     0 

lastic) 
. OOOE+00 
.OOOE+00 
.OOOE+00 
.OOOE+00 
. OOOE+00 
. OOOE+00 


2 .500E- 
.0000E+ 

.000E+ 
.000E+ 
.000E+ 
.OO0E+ 
.000E+ 
.000E+ 


3  .  614E  +  03 
0  0 

lastic) 
.OOOE+00 
.OOOE+OO 
.OOOE+00 
.OOOE+00 
.OOOE+00 
.OOOE+00 


7.227E+ 
. 0000E+ 

.000E+ 
'.  000E+ 
.000E+ 
.000E+ 
.000E+ 
.000E+ 


OOOE+00 
0     .OOOOE+00    .0000E+00  2 

OOOE+00  .OOOE+00  .OOOE+00 

OOOE+00  .OOOE+00  .OOOE+OO 

OOOE+00  .OOOE+00  .OOOE+00 

OOOE+OO  .OOOE+00  .OOOE+OO 

OOOE+OO  .OOOE+OO  .OOOE+00 


OOOE+00 


.OOOE+00       .OOOE+OO 


OOOE+00 
0     .0000E+00 


0000E+O0 


OOOE+00 
OOOE+00 
OOOE+00 
OOOE+00 
OOOE+00 
OOOE+00 


.OOOE+OO 
.OOOE+00 
.OOOE+OO 
.OOOE+00 
.OOOE+OO 
.OOOE+00 


.OOOE+00 
.OOOE+OO 
.OOOE+OO 
.OOOE+00 
.OOOE+00 
.OOOE+00 


.OOOE+00 

0     .0000E+00    .0000E+OO  1 

.OOOE+00  .OOOE+00  .OOOE+00 

.OOOE+00  .OOOE+OO  .OOOE+00 

.OOOE+00  .OOOE+00  .OOOE+00 

.OOOE+00  .OOOE+00  .OOOE+00 

.OOOE+OO  .OOOE+00  .OOOE+OO 

.OOOE+00  .OOOE+OO  .OOOE+00 


03    3.970E+02       .OOOE+00 

00  0     .OOOOE+00    .OOOOE+00 


00 


OOOE+OO 


00  .OOOE+00 

00  .OOOE+00 

00  .OOOE+OO 

00  .OOOE+00 

00  .OOOE+00 


OOOE+OO 
OOOE+00 
OOOE+00 
OOOE+00 
OOOE+OO 
OOOE+00 


000E+0C 
OOOE+00 
OOOE+00 
OOOE+00 
OOOE+00 
OOOE+00 


.OOOE+OC 
.  OOOE+00 
.  OOOE+00 
.  OOOE+00 
.OOOE+00 
.OOOE+00 


6 .275E+ 
OOOOOOO 
3291358 
0323486 
4315338 
3484497 
S879516 
OOOOOOO 
3315200 
0418853 
4219970 
3484497 
5879516 


04  1.255E+ 


OE+OO 
9E+01 
3E+01 
1E+01 
1E+01 
6E+02 
OE+OO 
8E+01 
8E+01 
7E+01 
1E+01 
6E  +  02 


.12 
.12 
.12 
.12 
.12 
.12 
.12 
.12 
.12 
.12 
.12 
.12 


05  8.880E+02 

65625000000E+ 

65625000000E+ 

65625000000E+ 

65625000000E+ 

65625000000E+ 

65625000000E+ 

09375000000E 

09375000000E+ 

09375000000E+ 

09375000000E+ 

09375000000E+ 

09375000000E+ 


.OOOE+00  .OOOE+OO  .OOOE+OO 

.OOOE+00  .OOOE+00  .OOOE+00 

.OOOE+OO  .OOOE+OO  .OOOE+OO 

.OOOE+00  .OOOE+00  .OOOE+00 

.OOOE+OO  .OOOE+00  .OOOE+00 

.OOOE+00  .OOOE+00  .OOOE+00 


.OOOE+00 

03  -.1681300163269E+02  5. 

03  -  .1660600662231E+02  0. 

03  -  .1599011611938E+02  0. 

03  -.1498049640G56E+02  0. 

03  -.1360200595856E+02  0. 

03  -.1188858795166E+02  0. 

03  -  .1681300163269E+02  5. 

03  -.1660600852966E+02  0. 

03  -.1599011707306E+02  0. 

03  -.1498049640656E+02  0. 

03  -.1360200691223E+02  0. 

03  -.1188858795166E+02  0. 
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913  0.  .12S9073257446E+02 

914  0.  .  1386673641205E+02 

915  0.  .1480129528046E+02 

916  0.  .1537139701843E+02 

917  0.  .1556300067902E+02 

918  0.  . 1537139511108E+02 

919  0.  .1480129432678E+02 

920  0.  .1386673641205E+02 

921  0.  . 1259073352814E+02 

922  1.  .0OOOO00O0O000E+00 

923  0.  .2434S8962440SE+01 

924  0.  .4809231758118E+01 

925  0.  .7065454959869E+01 

926  0.  .9147703170776E+01 

927  0.  .1100470352173E+02 

928  1.  .0OOOO0OO0O0OOE+OO 

929  0.  .2434589624405E+01 

930  0.  .4809231758118E+01 

931  0.  .7065454959869E+01 

932  0.  .9147703170776E+01 

933  0.  . 1100470352173E+02 

934  0.  . 1259073257446E+02 

935  0.  . 1386673641205E+02 

936  0.  . 1480129528046E+02 

937  0.  . 1537139701843E+02 

938  0.  . 1556300067902E+02 

939  0.  .1537139S11108E+02 

940  0.  . 1480129432678E+02 

941  0.  .1386673641205E+02 

942  0.  . 12S9073352814E+02 


.1153125000000E+03 
.1153125000000E+03 
-1153125000000E+03 
.1153125000000E+03 
.1153125000000E+03 
•11S3125000000E+03 
.1153125000000E+03 
.115312S000000E+03 
.1153125000000E+03 
.1209375000000E+03 
.1209375000000E+03 
.  120937S000OOOE+03 
.1209375000000E+03 
.1209375000000E+03 
.1209375000000E+03 
.1209375000000E+03 
.1209375000000E+03 
.1209375000000E+03 
.1209375000000E+03 
.  1209375000000E  +  03 
.1209375000000E+03 
.1209375000000E+03 
.1209375000000E+03 
.1209375000000E+03 
.  1209375000000E+03 
.1209375000000E+03 
.1209375000000E+03 
.  1209375000000E+03 
.  1209375000000E  +  03 
.1209375000000E+03 


.9147703170776E+01  0. 

.706S455436707E+01  0. 

.4809232711792E+O1  0. 

.2434590816498E+01  0. 

.0000000O00000E+0O  0. 

.2434589624405E+01  0. 

.48092317S8118E+01  0. 

.70654549S9869E+O1  0. 

.9147703170776E+01  0. 

.1556300067902E+02  5. 

.1537139511108E+02  0. 

.1480129432678E+O2  0. 

.138667364120SE+02  0. 

.1259073352814E+02  0. 

.11004703S2173E+02  0. 

.1556300067902E+02  5. 

.1537139511108E+02  0. 

.  1480129432678E  +  02  0. 

.  1386673641205E  +  02  0. 

.  1259073352814E  +  02  0. 

.  1100470352173E  +  02  0. 

.9147703170776E+01  0. 

.7065455436707E+01  0. 

.4809232711792E+01  0. 

.2434590816498E+01  0 

.0000O00000OOOE+OO  0. 

.2434589624405E+O1  C 

.4809231758118E+01  0 

.7065454959869E+01  0 

.9147703170776E+01  0 


ELEMENT  CARDS  FOR  BEAM  ELEMENTS 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 


553 
554 
555 
556 
557 
558 
559 
560 
561 
562 
563 
564 
565 
566 
567 
568 
569 
570 
571 
572 
573 
574 
575 
576 
577 
578 
579 


S54 

552 

555 

552 

556 

552 

557 

552 

558 

552 

559 

552 

560 

552 

561 

552 

562 

552 

563 

552 

564 

552 

565 

552 

566 

552 

567 

552 

568 

552 

569 

552 

570 

552 

571 

552 

572 

552 

573 

552 

574 

552 

575 

5S2 

576 

552 

577 

552 

578 

552 

579 

552 

580 

552 
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2.500E-01  2.500E-01    2.500E-01  2.S0OE-01 

839  7              920              921  383  364 
2.5O0E-O1  2.50OE-01    2.S0OE-O1  2.5OOE-01 

840  7               921              912  132  383 
2.500E-01  2.500E-01    2.500E-01  2.500E-01 

841  7                    7              922  923  8 
2.500E-01  2.500E-01    2.500E-01  2.500E-01 

842  7                    8              923  924  9 
2.500E-01  2.500E-01    2.5O0E-O1  2.500E-01 

843  7                    9              924  925  10 
2.500E-01  2.500E-01    2.500E-01  2.500E-01 

844  7                 10              925  926  11 
2.500E-01  2.500E-01    2.500E-01  2.500E-01 

845  7                 11               926  927  12 
2.500E-01  2.500E-01    2.500E-01  2.500E-01 

846  7              928              121  122  929 
2.500E-01  2.500E-01    2.500E-O1  2.500E-01 

847  7               929              122  123  930 
2.500E-01  2.500E-01    2.500E-01  2.500E-01 

848  7               930              123  124  931 
2.500E-01  2.500E-01    2.500E-01  2.500E-01 

849  7               931              124  125  932 
2.500E-01  2.500E-01    2.500E-01  2.500E-01 

850  7               932               125  126  933 
2.500E-01  2.500E-01    2.500E-01  2.500E-01 

851  7              927              934  230  12 
2.500E-01  2.500E-01    2.500E-01  2.500E-01 

852  7               934               935  249  23C 
2.500E-01  2.500E-01    2.500E-01  2.500E-01 

853  7               935              936  268  249 
2.500E-01  2.500E-01    2.500E-01  2.500E-01 

854  7               936               937  287  266 
2.500E-01  2.500E-01    2.500E-01  2.500E-01 

855  7               937              938  306  287 
2.500E-01  2.500E-01    2.500E-01  2.500E-01 

856  7               938              939  325  306 
2.500E-01  2.500E-01    2.500E-01  2.500E-01 

857  7               939              940  344  325 
2.500E-01  2.500E-01    2.500E-01  2.500E-01 

858  7               940              941  363  344 
2.500E-01  2.500E-01    2.500E-01  2.500E-01 

859  7               941              942  382  363 
2.500E-01  2.500E-01    2.500E-01  2.500E-01 

860  7               942              933  126  382 
2.500E-01  2.500E-01    2.500E-01  2.500E-01 

* 

* LOAD    CURVE    DEFINITIONS 

* 

1  2 
.000E+00       .OOOE+OO 

1.000E+00       .OOOE+00 

2  2 
.00OE+OO  2.000E-06 

1. OOOE+OO  2.000E-06 
* 

PRESSURE    BOUNDARY   CONDITION   CARDS 


1 

1 

7 

8 

1 

7 

13 

14 

1 

13 

19 

20 

2-1. OOOE+00 -1 . OOOE+00 -1  .  OOOE  +  00- 1 . OOOE+OO 

8-1. OOOE+00 -1 . OOOE  +  00 -1  .  OOOE  +  00 -1 . OOOE+00 

14-1. 000E+ 00-1. OOOE  +  00 -1  .  OOOE  +  OO- 1 . OOOE  +  0  0 
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FLUMAS    DATA   FOR    BEAMSHELL   MODEL 
beam. flu   beam.geo    strnam.dat    beam 
T   T   F  T 


TFT 

F 

F   F   T 

T 

F    F    F 

F 

F   T    F 

F 

F    F    F 

F 

F    F    F 

F 

DYNA 

0    942 

0 

450 

10    0 

Oil 

2 

0 

33 

0 

9.356E-05       57600.0 

5 

1.       1. 

0  .       -1. 

0.       0. 

3 

0.0 

16.813 

180 

1 

16 .812 

\ 

0.0            181 

360 

1 

0.0 

0.0            431 

520 

1 

0 

553    585 

16.813       16 

.813 

2 

553    553 

0.0            16 

.813 

1 

553    585 

16.813       16 

.813 

32 

daa 
$ 

$ 
$ 
$ 
$ 
$ 
$ 
S 
$ 
$ 
$ 
$ 

s 
s 
$ 
s 
s 
$ 
s 


FLUNAM      GEONAM      GRDNAM       DAANAM 


PRTGMT 

EIGMAF 

PCHCDS 

FRWTFL 

RENUMB 

PRTCOE 

OCTMOD 

MAINKY 

NSTRC 

NBRA 

NHAS 

NSEG 

NCIR 

RHO 

NVEC 

CQU) 

DHALF 

NSRADI 

RAD1 


PRTTRN 
TWODIM 
NASTAM 
FRWTGE 
STOGMT 
STRMAS 
CAVFLU 

NSTRF 

NCYL 

NHAF 


CEE 


CXHF 
RAD2 


PRTAMF 
HAFMOD 
STOMAS 
FRWTGR 
ROTGEO 
SPHERE 
FRWTFV 

NGEN 
NCAV 
NHAI 


CYHF 
JBEG 


CALCAM 
QUAMOD 
STOINV 
FRESUR 
ROTQUA 
ROTSYM 
INTCAV 

NGENF 

NFUN 


CZHF 
JEND 


NSORDER 
Nl       N2 


Rl 


R2 


JINC 


NSET 


AUGMAT    DATA    FOR 
strnam.dat    beam. 
F   F   F   F 

BEAMSHELL 
flu    beam. 

. MODEL 
geo   beam 

.  pre 
$ 

$       STRNAM      FLUNAM      GEONAM       PRENAM 
FRWTGE       FRWTST      FRWTFL       PLNWAV 

F   F    F   T 

$ 

FLUSKY 

DAAFRM 

SYMCON 

DOFTAB 

F    F    F  T 

$ 

PRTGMT 

PRTTRN 

PRTSTF 

PRTAUG 

F   F    F    F 

$ 

MODTRN 

STRLCL 

INTWAT 

CFADYN 

F 

$ 

LUMPFM 

DYNA 
0.50 
942    2826    3 

1 

3 

$ 
$ 
$ 
$ 

MAINKY 
DAA2 
NSTR 
NSETLC 

NSFR 

NFRE 

NFTR 

0    1    450    1 

$ 

NDICOS 

JSTART 

JSTOP 

JINC 

3 

$ 

NUMCON 

0       1      86 

5 

$ 

ICON 

NSTART 

NSTOP 

NINC 

0      91    176 

5 

0    400    472 

6 

105 


r=d3dump 

TIME  INTEGRATION  DATA 

beam. pre  beam.pos 


FOR  MODEL  BEAMSHELL 
$   PRENAM   POSNAM 


beam. rst 

$ 

RESNAM 

F  F  F 

$ 

BUBPUL 

REFSEC 

FLUMEN 

T  F  F  F 

$ 

EXPWAV 

SPLINE 

VARLIN 

PACKET 

F  T  F  F 

$ 

HYPERB 

EXPLOS 

DOUBDC 

VELIMP 

1 

$ 

NCHARG 

0.0 

$ 

HYDPRE 

0.0  -25. 3125  326  .4 

$ 

XC    YC 

ZC 

0.0  -25.3125   16.813 

$ 

SX    SY 

.SZ 

101 

$ 

JPHIST 

1.0 

$ 

PNORM 

14 .0E-06 

$ 

DTHIST 

1 

$ 

CHGTYP 

60.0   25. 8   10.0 

$ 

WEIGHT 

SLANT 

CHGDEP 

1  0 

$ 

NT  I  NT 

NCHGAL 

0.0   2.0000E-06 

$ 

STRTIM 

DELTIM 

5000   5000 

S 

NSAVER 

NRESET 

0  0  0  0 

s 

LOCBEG 

LOGRES 

LOCWRT 

NSTART 

F  F 

s 

FORWRT 

STBDA2 

F 

$ 

DISPLA 
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